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Chapter 1 
Introduction 


Magnetic resonance imaging (MRI) is a powerful non-invasive diagnostic tool in 
medicine because of its non-radioactive nature and high tissue contrast compared to 
other medical imaging modalities like X-ray CT. Magnetic resonance (MR) is based on 
the interaction between a system of atomic nuclei and an external magnetic field. A 
typical MR experiment consists of perturbing this system from its equilibrium 
configuration, and observing how it relaxes back to the equilibrium. During this 
relaxation process an induced electric signal is measured. This signal depends on three 
quantities proton spin density (PD or p), the interactions of the nuclei with their 
surrounding molecular environment (T i relaxation time), and those between close nuclei 
(T 2 relaxation time). Based on a suitable set of measured electric signals, a function of 
p, Ti, and T 2 is computed which is called MR image. In this thesis, the word 'tissue 
parameters' is used to represent all proton density, Ti and T 2 relaxation times. 

In nuclear magnetic resonance (NMR), the estimation of relaxation times dates 
back to shortly after discovery of the NMR phenomenon itself. These relaxation times 
estimation has been used extensively in NMR spectroscopy largely for the determination 
of the longitudinal relaxation time, Tj. Since the advent of magnetic resonance imaging 
in the early 1980s, parameters estimation has been extended to in vivo applications. 
NMR parameters have been proven useful in a variety of applications, including 
optimization of NMR pulse sequences, understanding of the contrast mechanisms in 
MRI, and identification of normal and disease processes. 

In biology, much of the interest in tissue relaxation times has started from the 
observation by R. Damadian (1971) that spin-lattice relaxation times are elevated in 
tumorous tissues (see P.G. Morris (1986)). This was demonstrated using two malignant 
rat tumors, Novikoff hepatoma and Walker sarcoma, which had Ti values at 24 MHz of 
736 and 826 ms respectively, compared with values in normal tissues extending over the 
range of 257 to 595 ms. Confirmation that this was a general effect followed rapidly 
with measurements of Ti values in animal and human tissues. Soon after Damadian’ s 
initial discovery, C.F. Hazlewood, D.C. Chang, D. Medina, G. Cleveland and B.L. 



Nichols (1972) showed that it was possible to detect cancer at early stage of 
development. They were able to distinguish normal, preneoplastic and tumourous 
murine glands not only on the basis of elevated Ti values but also by way of their 
increased T 2 values and diffusion constants. 

The specific goal of tissue characterization using an imaging technique is the 
possibility of discrimination between a lesion and adjacent normal tissues, between 
benign and malignant neoplasias, and between primary tumors and metastases using 
special analysis procedures. The objective is to improve diagnosis and prognosis and 
optimize therapy planning. Intrinsic multi parameter dependence is a peculiarity of MRI, 
which provides both high tissue contrast flexibility and a possible special source of in 
vivo tissue characterization. A primary role in the expression of these features is played 
by MR spin-lattice and spin-spin relaxation times of tissues in the body (see, F. Podo 
(1990)). An objective evaluation of the clinical significance of tissue characterization by 
MRI can, however, only derive from quantitative measurement of the physical 
parameters involved, together with their correlations with relevant clinical, biochemical, 
and histopathological data. 

This introductory chapter describes some basic concepts of MRI, two important 
pulse sequences viz., saturation recovery spin echo (SRSE), and inversion recovery spin 
echo (IRSE) and is followed by a summary of contents of the rest of the chapters of the 
thesis. 


1.1. Basics of Magnetic Resonance Imaging 

For the basics of Nuclear Magnetic Resonance (NMR) imaging, a number of excellent 
texts, e.g., T. Farrar and E. Becker (1971), P. T. Callaghan (1991), P. Mansfield and P. 
Morris (1982), C.L. Partain, R.R. Price, J.A. Patton, M.V. Kulkarni and E.V. James Jr. 
(1988), and D.D. Stark and W.G. Bradley Jr. (1999) are available. 

A simplified block diagram of a computer-controlled MR system is shown in the 
Figure 1.1. The computer and accessories constitute the command and control center of 
the MR system, performing shaping and timing the radiofrequency pulses, turning the 
gradients on and off, selecting region of interest (ROI) of the scanning object, 
instructing the receiver to collect data, manipulation, storage, retrieval, multi-format 
processing and providing the data for diagnosis. 




Figure 1.1: Simplified block diagram of a representative MR system 


1.1.1. Precession 

In presence of an applied magnetic field Ho, the magnetic momentum vectors 
experience a torque and precess about the axis of the magnetic field with a rate given by 
the Larmor's relationship 

=(cy/2;r) = (y/2;r)Ho, (1.1.1) 

where ft is the resonance frequency, (O is the angular frequency and y is the 
gyromagnetic ratio. The gyromagnetic ratio is equal to the ratio of the magnitude of the 
magnetic moment of the nucleus to the magnitude of its spin angular momentum. 









Gyromagnetic ratios, relative detection sensitivities, and isotopic abundances for 
a few diagnostically relevant nuclei are listed in Table 1.1 (M.L. Wood and F.W. Wehrli 
(1999)). The proton, which is the major isotope of hydrogen, is favored for MRI because 
of its abundance and sensitivity. 


Nucleus 

Relevant 

Relative 

Gyromagnetic 


Abundance (%) 

Sensitivity 

Ratio (MH,/T) 

‘H 

99.98 

1 

42.58 


0.015 


6.53 

‘^C 

1.11 

0.016 

10.71 

'^F 

100 

0.830 

40.05 

^■'Na 

100 

0.093 

11.26 

31p 

100 

6.6 X 0.066 

17.23 


93.1 

5.08 X 0.000508 

1.99 


Table 1.1: Magnetic resonance properties of some diagnostically relevant nuclei 


1.1.2. Radiofrequency (RF) Field 

MR can be detected only if transverse magnetization (magnetization perpendicular to 
Ho) is created because this transverse magnetization is time dependent and thus, 
according to Faraday's law of induction, can induce a voltage in a receiver coil. The 
longitudinal magnetization in thermal equilibrium is static and therefore does not meet 
the criteria for magnetic induction. Transverse magnetization is created when an RF 
field of amplitude Hi, rotating synchronously with the precessing spins, is applied 
(Figure 1.2(a)). When this RF field is directed perpendicular to the main field, the net 
magnetization rotates away from its equilibrium orientation. If Hi field rotates the net 
magnetization by 90°, all of the longitudinal magnetization is converted to transverse 
magnetization, as in Figure 1.2(b) (M.L, Wood and F.W. Wherli (1999)). 

1.1.3. Free-Induction Decay (FID) 

Once the RF pulse is removed, the magnetization processes about the static magnetic 
field at the Larmor frequency. The precessing magnetization can be detected as a time- 
varying electrical voltage across the ends of coil of wire placed around the object. The 




























(a) 


(b) 


Figure 1.2: Tipping of longitudinal magnetization into transverse magnetization by a 
magnetic field Hi, associated with an RF pulse, (a) For appreciable tipping to occur, the Hi 
field must rotate synchronously with the processing transverse magnetization. This condition 
is met by circularly polarized RF pulse at the Larmor frequency, (b) All of the longitudinal 
magnetization is tipped into the transverse plane by an RF pulse that is strong enough and 
applied long enough to tip the magnetization 90° (90-degree RF pulse) 


magnetization also decays exponentially with a time constant T 2 . A simple model for 
this induced voltage is (M.L. Wood and F.W. Wherli (1999)) 

V = K Mq txpiilTtfj) exp(-r /T^), 

where x: is a constant , i = (-1)'^^, and Mo is magnetic momentum at thermal 
equilibrium. The induced voltage has the characteristics of a damped cosine and hence 
is also called free-induction decay (FID). 

1 .1 .4. Proton Density (PD or p) 

The overall strength of the signal is decided by the proton density (PD or p) that makes 
up the macroscopic magnetization. The proton density relevant in MRI is the number of 
mobile protons per unit volume in contrast to the total number of protons (that include 
stationary protons). Typically, for example, protons that belong to the organic matter in 
dense bone do not contribute to the signal that comprises the MR imaging experiment. 
They produce an MR signal, except that its characteristics are such that it is not relevant 
to the MR imaging process - more specifically, the signal from protons belonging to 
solids decays at too fast a rate (~ tens of microseconds) compared to mobile protons 



(hundreds of milliseconds) for any meaningful data to be collected. For this reason, 
solids, with an absence of mobile protons appear dark in an MR image. However, the 
lack of signal does not render such information useless - a lot can be surmised from the 
appearance of the anatomy surrounding the solid structure. Typical structures that 
appear dark due to a lack of protons are cortical bone, sinus passages, lung parenchyma 
and calcium deposits. Changes in proton density occur from changes in the cellular 
composition in tissues that accompanies disease and will be reflected in the signal 
intensity (S. Sinha M.A. Thomas, A.S. Huda, and R.B. Lufkin (2001)). 

1.1.5. Spin-Lattice Relaxation Time (Ti) 

RF stimulation causes nuclei to absorb energy, lifting them to the excited state. The 
nuclei can return to the ground state only by dissipating their excess energy to their 
surroundings, which is called the lattice. The process, which is known as the spin-lattice 
relaxation or Ti relaxation describes the recovery of the longitudinal magnetization 
toward its equilibrium value. The relaxation of the excited spin requires an interaction 
with a time-dependent perturbation, typically a dipole-dipole interaction with a 
fluctuating electromagnetic field at the Larmor frequency. Such perturbations or 
fluctuating fields are produced by rotation and translation of the nuclei in the molecules 
of the lattice undergoing Brownian motion. The average frequency of this Brownian 
motion depends on the size of the molecules in the lattice. Small molecules re-orient 
more rapidly than larger molecules. The frequency of rotation in medium size 
molecules, such as lipids, is closest to typical Larmor frequencies. Therefore the 
magnetization associated with lipids relaxes faster than that associated with pure water 
or much larger molecules, such as proteins. Moreover, Ti relaxation times depend on the 
magnetic field strength because the latter affects the Larmor frequency. 

1.1.6. Spin-Spin Relaxation Time (T 2 ) 

Before the RF pulse is switched on, the individual magnetic spins (protons) are on the 
average aligned with the main magnetic field Ho, and in random distribution of phase 
with respect to each other. The RF pulse not only tips the magnetization away from the 
z-axis but also puts them in phase coherence with each other. This coherent 
magnetization tipped into the transverse or xy plane, precesses about the Ho (z-axis), 
inducing the MR signal at the Larmor processional frequency. The transverse 



magnetization decays because its component magnetic moments get out of phase as a 
result of their mutual interaction and thereby are exposed to static as well as slowly 
varying local magnetic field fluctuations. These fluctuations cause a spread of the 
magnetic field around the central value of Ho, causing each to be exposed to a slightly 
different magnetic field, which in turn causes them to process at slightly different 
angular frequencies, some faster than /£, others slower. The net effect is that the ordered 
phase coherent state slowly disappears, with the spins fanning out in the xy plane. This 
causes the induced MR signal to decay at a rate governed by a factor called T 2 or, the 
spin-spin relaxation time since it is caused primarily by intrinsic fluctuations in 
magnetic field arising from interactions between adjacent spins, or alternatively, the 
transverse relaxation rate, since it takes place in the transverse plane. 

Unlike Ti relaxation, no energy is transferred from nuclei to the lattice in T 2 
relaxation. Nuclei in the excited and ground state may exchange energy with each other. 
However, in biological tissues the main contribution to T 2 relaxation is from the 
relatively static magnetic field from neighboring protons. Large molecules such as those 
of protein, which tend to reorient more slowly than small molecules, promote T 2 
relaxation and have shorter T 2 times. In contrast, smaller molecules such as water, at 
least at body temperature, have very rapid Brownian motion that average out any field 
fluctuations, or produce high frequency components. The result is that water has very 
large T 2 values, which implies that once a transverse spin-spin coherence is built up in 
an ensemble of spins, this coherence takes a long time to dephase - water and fluids 
such as cerebrospinal fluid (CSF) and edema have long T 2 ’s. The magnetic field 
strength influences T 2 much less than Ti, at least under the conditions encountered in 
MRI. 

Table 1.2 shows the Ti and T 2 relaxation times and proton spin density values of 
some brain tissues measured at 1.5 Tesla (R.E. Hendrick (1999)). Here, the proton spin 
density values are normalized with respect to the CSF proton density value. 


Tissue 



P 

White matter 


67 


Gray matter 

760 

77 


Edema 

900 

126 

0.86 

Cerebrospinal fluid 


280 

1.00 

Upper cystic fluid (extracellular methomoglobin) 

1080 

215 

1.12 

Lower cystic fluid (intracellular deoxyhemoglobin) 

720 

43 

1.02 

Methomoglobin 



0.96 


Table 1.2: Brain tissue parameters measured in selected section at 1.5 T 

























1.1.7. Bloch Equations 


F. Bloch (1946) found that the motion of the macroscopic magnetization in the presence 
of an applied magnetic field could be explained in terms of phenomenological 
differential equations. The classical equation of motion of a magnetic moment in a 
magnetic field H states that the rate of change of angular momentum of the spinning 
nucleus, p, depends upon the torque, q x H, exerted on the magnetic moment by the 
applied field and can be represented as: 

Jp 


dt 


= /zxH. 


Multiplying both sides of this equation by y and using )i=Yp, we obtain 

djj, dp „ 

-p = r-f = 7;^xH. 

dt dt 

If M is the vector sum of the p, 's, then by summing equation vectorially over all p we 
obtain the following relation for macroscopic magnetization: 

cM 


dt 


■ = 7MxH. 


( 1 . 1 . 2 ) 


During the RF pulse, the magnetic field is 

H(r) = //|(r) cos cun' + Hi{t)sin(Ot j + Hk . 

Here, H = Hq + H, where h is the contribution from the gradients. Thus the equation 

(1.1.2) can be written component wise as: 

dM. 


dt 

dM 

dt 

dM 

dt 


- = YM yli -y M ,H^{t)sinCL)t , 

- = YM coscot -Y M , 

= 7 sin cot -Y coscot . 


Let M, =ucoscot -vsincot and M„ =usincot + vcoscot . 

X y 

Then the differential equation satisfied by u and v is 




^ u 


■ 0 

Y H +CU 

0 ^ 

where X = 

V 

and A = 

-(yH +co) 

0 

YHx 


M, , 
^ J 


[ 0 

-yh, 

0 

y 


Therefore X{t) = exp(Ar)X(0) . 



Now, 


-A 0 

|A-AI| = 0=> -(y/Z+oj) -A y/Z; =0=> A^ +(y^//- +(y// +a))-)A = 0. 
0 -yZZ, -A 

Therefore the three eigen values of the matrix A are 
A, =0;A2,3 =± +{yH +0}f . 

( 0 a 0^ 

If we write A= -a 0 b where a=yH +0) and then 

0 -b 0 


0 1 and AU* 

T— 0 

la-+b~ 


0 

a- +b- 
0 


"0 0 

Therefore the transformed matrix is 0 0 

0 


-^|a^ +b~ 


Write C 


c 0 0 1 


= -c^ 


0 -1 

1 0 


;C'=c 


and so on, 


where c = sa‘' +b^ . 


Therefore exp(a) = |;[^.^ 

2n\ (2n + l)! 


£^ 2 / 1 ! 10 1) (2n + l)! 


0 -1 

1 0 


f\ O'] (0 -1 

cos ct + sin ct 

0 1 10 


cos ct - sm ct 
sinct cos ct 


Therefore X(0 = exp(At)X(0) 

= exp(UDU*)X(0) 


^10 0 - 
U 0 cos ct -sincf U*X(0). 
0 sinct cosct 



Case (I): If the initial condition is X(0) = (O 0 M^J , then the solution will be 


^ u ^ 




abMf^ 

0 TT 
a" ->rb~ 


(l-cosct) 


bM^ 


0 


i 


-sinct 


a^+b^ 

2.72 


^ (q- +&-^coscO 
a“ 


Therefore 


•^v(0 = - - ■ Mq (1 - cos ct) cos (Ot — ;■ ■ - Mq sinct sin cut , 

M ^{t) = Mq (1 - cos ct) sin (tit + ■ i -~—- Mq sin ct sincot , 


M(t) 


a^+b- 


^|a^ +b^ 


(a~ +b^cosct ) , 


■\ja' +b~ 


Case (11): If the initial condition is X(0) = (m^.(0) M^(0) M,(0))^ , then 


^ u ^ 


\^z J 


^ j{b~ + a' cosct)M ^{0) + - Y^-^--=^ sinct M y{0) + ^ ^ , (l-coscr)M,(Q) 

■\]a'^+b^ a" -^b- 


0 , n 

a~ + b 




^ 7 sin ct M (0) + cos cr M ^ (0) + - ^-L=.=^ sin ct M , (0) 




\a' +b^ 

(l-cosct)M^(O) — , - •- sinctMy(Q) + -y-^ — +b' cos ct)M^iO) 


y 


" 7 2 

a~ +b 


i 


+b'^ 


a^+b 


Therefore 

f yi ^2 ^ 

r COS tyt + -T cos ct cos cot + , sin ct sin cot 

ct'+Z)- a +& yla^+b^ 

+ M (0)\ sin cot coscot - cos ct sin cot 

+ M , (0) i ■ cos cot — cos ct cos cot - ■ sin ct sin cot 

^ a' + 6' ^a^+b^ 

\ b^ CL^ ■ a . 

M(t) = M(0)\ — 7 sin cot + — -Y cos ct sin cot — /--- -- sin ct cos cot 

" a^+b^ sla^+b^ 

+ M (0)|— p=^=sinctsincot-cosctcoscot 

" [47^ 

+ M (0)- sin cot — coset sin cot + , ^ ■ sin ct cos cot 

a^+b^ a^+b^ 444^ 




M,(r) = M,.(0) 


ab 


a^+b^ 


(l-coscf)-M (0) 


+b‘ 


rsincr + M ,(0) 


+b^ 


(a^ +b" coscr). 


Since the above equations do not account for relaxation, F. Bloch et al. (1946) assumed 
that spin-lattice and spin-spin relaxations could be treated as first order processes with 
characteristic times Ti and T 2 respectively. Mx and My decay back to their equilibrium 
value of zero, while Mz returns to its equilibrium value of Mq. 

In the absence of RF and considering Ti , T 2 relaxation, Bloch equations become 


dM^ 

dt 


M,. 


= YM^H- ^ 


dM^ 

dt 


-YM^H 



dM,_ M,-Mo 
dt r, 

The third equation of the above coupled system of equations can be solved 
independently and the solution with the initial condition M,(0) is 

= Mo + (M,(0)-Mo)exp(-r/r,) . 


First and second equations of the above system can be written in the matrix form as 

-^X(0 = AX, 
dt 


where X^ 




and A 


>2 

■yh 


Solution will be 




YH 

1 



X(0 = exp(Ar)X(0) 

f 

= exp(-r/Tj 

V 


cos Yfit 
- sin YHt 


sin YHt 
cosYHt 


MM'' 

M/0)^ 


Therefore, 

(r) = exp(-r / 72 ) {M^ (0) cos YHt + My (0) sin YHt), 


My{t) = exp(-r / T 2 ) {M y (0) cos YHt-M^ (0) sin YHt ] , 
M^ (r) = Mq + (M^ (0) - Mq )exp(-r IT{). 



In actual MR experiment, due to local field inhomogeneities, the transverse 
magnetization relaxes at a faster rate T2* instead of T2. Considering spatial coordinate x 
= (x,y,z), and T2*, the above solutions can be written as 

M,(x,r) = exp(-r/r2*){M_,(x,0)cosy//r + M^(x,0)siny/fr} , (1.1.3a) 

M y (x, r) = exp(-r /T 2 ){My (x, 0) cos yHt-M^ (x, 0) sin yHt } , (1.1.3b) 

M , (X, r) = M 0 (X) + (M , (X, 0) - M 0 (x))exp(-t /T^). ( 1 . 1 .3c) 

By placing a conductive coil around the sample, with axis lying on the xy-plane 
by Faraday-Neumann’s law, it will induce, an electric signal s(t) in the coil. This is the 
only measurable quantity in an MR experiment. The signal s(t) is usually called free 

induction decay (FID) signal; it is a two-dimensional vector and thus it can be 

represented as a complex function 

s{t)o^ s^.{t) + is y{t) . 

Due to hardware constraints the signal s(t) is measured not in the laboratory frame, but 
in a frame, called signal detection frame, rotating around k with frequency (Or very near 
to the Larmor frequency (Oo = yHo- Following G. Sebastiani and P. Barone (1991), 
sit)o: expiio)^ -Q)Q)txj(M'^(x,t) + iM'yix,t))dx 

V 

= expi(co^ -ft)o)?J|M'j^.j,(x,0) xexp((i0i(x,r) + ^''(x))-r/r2*(x))rfx , (1-1-4) 

V 

t 

where 0, (x,r) = -yj H^ (x,r)dT : , and the symbol ' refers to Zeeman's frame. 

0 

Subsections 1.1.8-1.1.12 are mainly followed from G. Sebastiani and P. Barone (1991). 

1 .1 .8. Selective Excitation 

To select a slice parallel to the xy-plane at height zq, a radio frequency pulse lying in xy- 
plane will be applied, but with a time varying H (sine-shaped), together with a linear 
gradient along z-axis in the static field k , which is known as a selective pulse. In case 
of resonance, we will have 

Jlix,t) = H X ( [cos (Oq r i] - [sin cOq F] j) + [H 0 -H g ^ (z - Zo )] k , 

Acot 

where 7 = 0<t <t(A(0, H and are constants); the g,(z-Zo)*^ is called 

slice selection gradient. The application of a sine shaped radio frequency pulse cannot 



be implemented exactly because of its infinite time extension, but the truncated version 
of it {0<t <tp) will be used, including ripples in frequency spectrum and, as a 

consequence, in the slice profile. This effect can be reduced by using a sine function 
multiplied by a damped function, for example exp{-b{ t - Vi tpf) (b is a positive 
constant), so that the slice profile is similar to a Gaussian density. After the application 
of a 90° selective pulse, FID signal can be approximated by 



t 

T2{x,y,Zo) 


\ 

dxdy 

J 


(1.1.5) 


where S is the xy-projection of the intersection between V and the plane z = zo- 


1.1.9. Imaging Techniques 

The basic idea for the encoding process was introduced by P.C. Lauterbur (1973), to 
obtain the first MR images. After slice selection by a 90° selective pulse, the signal 
acquisition is performed in presence of linear gradient along x-axis in the field, i.e. 
H(x,0 = (7^o iSx a constant), so that the Larmor frequency is equal to 

Y^Hq + g a) ’ which is constant on the straight lines x = and is a linear function of x . 
This technique is known zs frequency encoding. 


1.1.10. Fourier Zeugmatography 


Zeugmatography is the name given by Lauterbur to what has now become MRI taken 
from the Greek Zeugma (to join together); referring to the combining of the static Hq 
field and its gradient with the NMR signal from tissue nuclei to produce a graphic 
image. A. Kumar et al. (1975) have introduced a process called phase encoding, which 
consists of the application, after the slice selection and before the frequency encoding, 
of a linear gradient y in the field, i.e. H(x,r) = {Hq + gyy)k for a time ty. During the 


time ty, the re-phasing field gradient -g^(z-Zo)k is also applied, but we can assume 
that this happens before the encoding process. The gradients and g^k are called 
readout frequency gradient and phase encoding gradient respectively. 

r' + r. 


s{t') oc jMo(x,y,Zo)xexp 




T2(x,y,Zo) 


+ iYig^xt'+gyyty) 


dxdy 


( 1 . 2 . 6 ) 



or 


/ f ^ +^■ 

s{t) oc I p(A-,y,Zo)xexp — ^ + iY{g^xt'+g^yt^) dxdy (1.2.7) 

5 (^7’2 (j:,j,Zo) ■ ■ 

Let us consider ensemble of signals s'^ (f ) , obtained from the experiments with g ^ , and 

g^. fixed and positive, and 0<ty <7^ (T^is a positive constant). For each of them the 
acquisition is formed for 0 < r'< F . 


F(m,^) = J js{t',ty)txg>i{Q)t'+^ty)dt'dty 


j p{x,y,Zo) 


sinjo) -Yg,x)T Y sin(g - Yg^.y)T^ 
oy-Yg,x I ^-YgyV 


\dxdy (1.2.8) 


As T,T^, — > oo, the two sine functions appearing in the Fourier transforms will tend, as 


distributions, to <5 (ft) - ,.x) and <5(^-ygyy), respectively, so that F{co,^) will 
approach to p(ft)/7g^,,^/yg^, ,Z q). It is clear that can not increase indefinitely 

without becoming large with respect to the values of T* . This is the limitation of 
Fourier zeugmatography. In order to overcome this problem W.A. Edelstein, J.M.S. 
Hutchison, G. Johnson, and T. Redpath (1980) proposed a technique called spin-warp. 


1.1.11. Spin-warp 

Let us consider the ensemble of signals obtained by the experiments with 

|g ,, < G^(Gyis a constant), and all other parameters are fixed. For each g^ value, the 
acquisition is performed for 0<t'<T. As before, Fourier zeugmatography, for any 
values t', g^suchthat g^ < G^ , we define the function s{t',gy) and apply the 

two dimensional inverse Fourier transform with respect to fand g^ . If T and /j,are 
small with respect to (*) values, we have 

oo 

F(ft),77) = J js(t',gy)expi{(ot'+r]gy)dt'dgy 



sin(ft) - yg^x)r 
ft)-yg^x 


oc 


smi'n-Ytyy)Gy 

v-r^yy 


\dxdy 


(1.2.9) 



1.1.12. Spin-echo 


spin-echo is similar to spin-warp: at time TE/2 from the end of the 90° selective pulse, a 
180° non-selective pulse is performed for a negligible time. Before this pulse, the field 
{Hq + g + g^x)k is applied which is substituted, after the pulse, by the field 

{Hq +g ^.j)k . Therefore, immediately before and after 180° pulse, the ^y-projection of 
the magnetization density in Zeeman’s frame has a phase equal to -yig^y + g^x)TEI2 
and y(g^,y + g^x)r£'/ 2 respectively, and its modulus is Mq{x)&x^{-TE l2T*_{x)) . In 


the following when the field + g ,.Ac)k is applied, a signal called spin-echo is 
produced at t = TE is proportional to 


JpW 


exd 


V V 


TE . TE 
T,{x) 2 


dx . 


For t > TE, signal decays because of x-gradient, spin-spin interaction and field 
inhomogeneities. Shifting time origin at t = TE, we have 


s{t) = J p(x)exp(-r£'/7’2)exp 


^ ^ r TE 


T^ix) 


dxdy. (1.2.10) 


5 V . 

As for spin-warp, the inverse Fourier transform is performed to decode the signal. 


1.2. Pulse Sequences 

One of the strengths of MR imaging is the myriad of measurement techniques, known as 
pulse sequences and parameter modifications, available for use. Even though there are 
many pulse sequences available, every pulse sequence must accomplish two basic tasks. 
First, transverse magnetization must be created using one or more radio frequency 
pulses, and this magnetization must then be encoded using gradient pulses so that the 
spatial positions of the tissues can be determined from the image. Second, the desired 
contrast between tissues must be created based on the time of the application and 
duration of the various RF and gradient pulses. 


1 .2.1 . Pulse-Sequence Timing Diagram 

A useful method for conveying the important characteristics of a pulse sequence is a 
graphical representation called the pulse sequence timing diagram. In general, at 



minimum five lines are needed to completely describe any pulse sequence (see M.A. 
Brown and R.C. Semelka (1999) and J.P. Mugler III (1999)). One for RF transmitter, 
one for each gradient and the last one for data acquisition. The first line, labeled RF, 
illustrates the time of application and waveforms for the RF pulses, in this case 90° 
(excitation) and 180° (refocusing) pulses. The next three lines depict the time of 
application and the waveforms for gradients applied along three mutually orthogonal 
axes. These are labeled GS, GP, and GF, for the section select, phase-encoding, and 
readout gradients, respectively, although other labels, for example Gx, Gy, Gz, can be 
used to indicate the assignment of the gradients for a particular section orientation. For 
the two dimensional (2D) acquisition depicted in the Figure 1.2, section-select gradients 
are applied in synchrony with the RF pulses, thus spatially localizing their effects, or in 
other words selecting the section of interest. The other gradient waveforms, applied 
along the phase-encoding and readout axes, spatially encode the magnetization within 
this section along the remaining two directions. The gradient waveform on the phase- 
encoding axis, containing a series of closely spaced horizontal lines, is the graphical 
representation of a gradient table: strength of the gradient is stepped through a series of 
values as the basic pulse sequence timing is repeated to collect the data required to form 
the image. With each repetition, another line of data is collected. The fifth line of the 
diagram, shows when the MR signal occurs. Although some pulse sequences generate a 
number of distinct signals, typically only the signals of interest are depicted. 

1.2.2. Saturation Recovery Spin-Echo (SRSE) Sequence 

Saturation recovery spin echo sequence is a routinely used sequence for medical 
diagnosis. It is also simply called as spin echo (SE) sequence. Depending upon the way 
it is implemented, some variants of this sequence are also available. For example, fast 
spin echo (FSE) sequence is a variant of SE. Pulse sequence timing diagram of SE is 
shown in Figure 1.3. 

In saturation recovery spin echo sequence, the magnetization density at the 
beginning of the 90° pulse of the first unit is Mo(x)k, and at the end of the pulse it lies 
on the jcy-plane with the magnitude Mo(x). At the beginning of the 180° pulse of the 
same unit the component along k is equal to Mo(x)[l-exp(-TE/2Ti(x))]. At the end of 
this pulse the component along k of the magnetization density will be -Mo(x)[l-exp(- 
TE/2Ti(x))]. At the beginning of the 90° pulse of the second unit, a time TR is elapsed. 
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Figure 1.3: Timing diagram of spin echo sequence 


and, since TR is long compared to T 2 (x), it follows from solutions of Bloch equation 
that the transversal components of the magnetization density will be zero, while Mz will 
be equal to 

Mo (x)[l + exp(-77? / T, (x)) - 2exp(-(77? -TEI2)IT^ (x))] 
with Mj.(x,0) = -Mo(x)[l-exp(-7’£'/2r,(x))]. 

At the end of the pulse, the magnetization density lies on the xy-plane with magnitude 
equal to the above value Mz. This is true for all the other units too. The signal 
expression will be obtained from equation (1.2.10) replacing p(x, y, zq) by 

Pix, y, Zo )& + exp(-rR / Tj (x, y, Zq))-2 exp(-(rR -TEI2)IT^ (x, y, Z q ))] ■ 

If TR»TE, then this last expression becomes p[l-exp(-TR/T^)]-, for two tissues A 
and B with the same spin density and different Ti, the difference of their contribution to 
the signal for unit volume is proportional to (exp(-TR/TiA)-exp(-TRyTiB)). The suitable 
value of TR is the one that maximizes this function and, as a consequence, the contrast. 

The transverse and longitudinal magnetic moment components Mt and Ml 
respectively which account for the relaxation process can be written as (J.P. Jones 
(1988a)): 


M ^ (0 = M (1 - exp(-r ITy)) + M (0) exp(-r / ) 




M-p(t) = Mj (0) exp(-? /T 2 ) (true decay) 

Mj (t) = M.J. (0) cxpi-t/T^ ) (FID decay) 

where Ml(0) is the value of Ml immediately following an RF pulse, t is any time after 
the pulse, and similarly for Mj. Longitudinal and transverse components of magnetic 
moment vector at different stages of first cycle of pulse sequence timing diagram are 
shown in Table 1.3. Table 1.4 shows the values for second cycle and these values are 
same for all other cycles. 


Event 

Ml 

Mt 

Start 

M 

0 

90° 

0 

M 

TE/2 

M(l-exp(l-TE/2Ti)) 

Mexp(-TE/2T2*) 

180° 

-M(l-exp(l-TE/2Ti)) 

Mexp(-TE/2T2*) 

TE/2 

M(l-2exp(-TE/2Ti)+exp(-TE/Ti)) 

Mexp(-TE/T 2 ) 



measure 

End 

M( 1 -2exp(-(TR-TE/2)/Ti)+exp(-TRn’i)) 

Mexp(-TE/T 2 )exp(-TD/T 2 *) 


Table 1.3: Longitude and transverse magnetic moment values at each step of the first spin-echo 


of SRSE 


Event 

Ml 

Mt 

Start 

M( 1 -2exp(-(TR-TE/2)/Ti)+exp(-TR/T,)) 

»0 

90° 

0 

M(l-2exp(-(TR-TE/2)/T,)+exp(-TR/Ti)) 

TE/2 

M(l-exp(l-TE/2Ti)) 

M(l-2exp(-(TR-TE/2)/T,) 

+exp(-TR/Ti))exp(-TE/2T2*) 

180° 

- M(l-exp(l-exp(TE/2Ti)) 

M(l-2exp(-(TR-TE/2)/Ti) 

+exp(-TR/Ti))exp(-TE/2T2*) 

TE/2 

M(l-2exp(-TE/2Ti)+exp(-TE/Ti)) 

M(l-2exp(-(TR-TE/2)/Ti) 

+exp( TRyTi))exp(-TE/T 2 ) 



measure 

End 

M(l-2exp(-(TR-TE/2)/Ti)-i-exp(-TR/Ti)) 

M(l-2exp(-(TR-TE/2)/Ti) 

-t-exp(-TR/Ti))exp(-TE/T 2 )exp(-TDAr 2 *) 


Table 1.4: Longitude and transverse magnetic moment values at each step of the second spin- 
echo of SRSE 











































Combinations of long values of TR with respect to those of Ti, and of the values 
of TE comparable with those of T 2 , are useful to weigh the image by T 2 , while 
combinations of short values of TE with respect to those of T 2 , and values of TR 
comparable with those of Ti, give images weighted by Ti. If long values of TR are 
associated to short values of TE produce images influenced essentially by spin density. 
These effects of TE and TR are shown in Table 1.5. 


Contrast Source 

TE Setting 

TR Setting 

T, 

As short as possible, limited by band 

width and noise (<15 ms) 

Comparable to TiS of the two tissues 

(better towards the shorter Ti) 

To 

Comparable to T 2 S of the two tissues 

(better towards the longer T 2 ) 

Long compared with the Tjs of the 

two tissues (>2000 ms) 

P 

As short as possible (<15 ms) 

Long compared with the TiS of the 

two tissues (>2000 ms) 


Table 1.5: Effect of TE and TR on image contrast in spin-echo imaging 


Figure 1 .4 shows the images of brain with different forms of contrast produced by TE 
and TR vales in spin echo sequence. Figure 1 .4(a) is obtained with TE and TR as 20 ms 
and 2500 ms respectively to produce a proton density weighted image. Figure 1.4(b) is a 
T 2 weighted image with TE = 80 ms and TR = 2500 ms. On T 2 weighted image white 
matter appears as dark, gray matter as gray and CSF as very bright. In contrast to this, 
on the Ti weighted image obtained with TE = 15 ms and TR = 650 ms (Figure 1.4(c) ) 
whiter matter appears bright and CSF as dark. All these three weighted images produce 
different contrasts and reveal excellent anatomic details. 
















(C) 

Figure 1.4: Images of the human head with different forms of contrast: (a) a spin density- 
weighted image, (b) a Tn weighted image and (c) a Ti weighted image. These different 
acquisitions can be seen to create different contrasts between white matter, gray matter and 
cerebrospinal fluid. They all reveal excellent anatomic detail. 


1.2.3. Inversion Recovery Spin-Echo (IRSE) Sequence 

Inversion recovery spin echo (IRSE) pulse sequence is used to weight the spin density 
by Ti too. In the literature of MRI this sequence is simply known as inversion recovery 
(IR). It consists of performing the same sequence as above, with an additional 180° non- 
selective pulse at time TI before the 90° selective pulse (180° non-selective pulse - 90° 
selective pulse - phase and frequency encoding - 180° non - selective pulse - frequency 
encoding with acquisition). The pulse sequence timing diagram of IR sequence is shown 
in Figure 1.5. 






Figure 1.5: Inversion recovery pulse sequence timing diagram 


At the beginning of the first 180° and of the 90° pulses of the first unit the 
magnetization density lies along k, so that at the end of the 90° pulse it lies on the xy- 
plane. If TR-TI is long with respect to the values of T 2 ; then, at the beginning(the end) 
of the first 180° pulse of the next unit, the transversal components of the magnetization 
density will be zero and Mz will be 

- Mo (1 - exp(-(r/? - r/) /Tj ) . 

The value of Mz at the beginning of the 90° pulse of the same unit, will be obtained by 
(1.1.3c) with the conditions 

M , (X, 0) = -Mo (x)[l - tx^{-{TR -T1)IT,)] 

at t = TI. Therefore we have that the magnetization density at the end of the 90° of each 
unit following the first one, lies in the xy-plane, with magnitude 

Mo [1 - 2 exp(-r/ / Ti ) + exp(-T/? / )] 

The expression of the signal is obtained from (1.2.10) by substituting p(x, y, zq) with 

p[l - 2 exp(-77 / Ti ) + exp(-ri? / T, )] . 

If TR is long compared to the values of Ti , then this last expression becomes 

p[l - 2 exp(-T/ / Tj )] . 

We note that, for two tissues with the same spin density but different Ti the maximum 
difference of their contributions to the signal for unit volume is twice that obtained with 
SRSE. 




Following the notation for Ml and Mt as given in the previous subsection, 
longitudinal and transverse components of magnetic moment vector at different stages 
of first cycle of pulse sequence timing diagram for IR sequence are shown in Table 1.5. 
Table 1.6 shows the values for second cycle and these values are same for all other 
cycles. 


Event 

Ml 

Mx 

Start 

M 

0 

180° 

-M 

0 

TI 

M(l-2exp(TI/Ti)) 

0 

90° 

0 

M(l-2exp(Tl/Ti)) 

TE/2 

M(l-exp(l-TE/2Ti)) 

M(l-2exp(TI/Ti))exp(-TE/2T2‘ ) 

180° 

-M(l-exp(l-TE/2Ti)) 

M( l-2exp(TI/Ti))exp(-TE/2T2* ) 

TE/2 

M(l-2exp(-TE/2T,)+exp(-TE/T,)) 

M(l-2exp(TI/Ti))exp(-TE/T2 ) 



^ measure 

End 

M(l-2exp(-(TD+TE/2)/Ti) 

+exp(-(TD+TE)/T,)) 

M( l-2exp(TI/ri))exp(-TE/T2) 

X exp(-TD/T 2 *) 


Table 1.6: Longitude and transverse magnetic moment values at each step of the first spin-echo 


of IRSE 


Event 

Ml 

Mx 

Start 

M( 1 -2exp(-(TD-)-TE/2)/Ti) 
-t-exp(-(TD+TE)/Ti)) 

0 

180° 

-M( 1 -2exp(-(TD-i-TE/2)/Ti ) 
+exp(-(TD-t-TE)/Ti)) 

0 

TI 

M(l+2exp(-(TR-TE/2)/Ti) 

-exp(-TR/Ti)-2exp(-TI/Ti)) 

0 

90° 

0 

M(l+2exp(-(TR-TE/2)/T,)-exp(-TR/T,)-2exp(- 

TI/Ti)) 

TE/2 

M(l-exp(l-TE/2Ti)) 

M(l-t-2exp(-(TR-TE/2)/Ti)-exp(-TR/Ti)-2exp(- 
TlfT,)) exp(-TE/2T2*) 

180° 

-M(l-exp(l-TE/2Ti)) 

M(l-i-2exp(-(TR-TE/2)/T0-exp(-TR/Ti)-2exp(- 
TI/Ti)) exp(-TE/2T2'') 

TE/2 

M(l-2exp(-TE/2T,) 

+exp(-TE/ri)) 

M(l+2exp(-(TR-TE/2)/Ti)-exp(-TR/T,)-2exp(- 
TITT,)) exp(-TE/T 2 ) 



4- measure 

End 

M(l-2exp(-(TD+TE/2)/r0 

-i-exp(-(TD+TE)/Ti)) 

M(l+2exp(-(TR-TE/2)/Ti)-exp(-TR/T,)-2exp(- 
TFTi)) exp(-TE/T 2 )exp(-TD/T 2 *) 


Table 1.7: Longitude and transverse magnetic moment values at each step of the second spin- 


echo of IRSE 

Inversion recovery accomplishes selective nulling of a specific tissue based on its T i 
rp.layatinn time and the chosen inversion time TI. Depending upon the choice of TI and 




the way the sequence is implemented, there are many variants of IR sequences routinely 
used in clinical diagnosis to achieve different contrast images . A few of them are short 
TI inversion recovery (STIR), fluid-attenuation inversion recovery (FLAIR), balanced 
inversion recovery (BIR), Fast inversion recovery (FIR), and magnetization prepared 
rapid acquisition gradient echo inversion (MP-RAGE). Following G.M. Bydder, J.V. 
Hajnal, and I.R. Young (1999), the effect of different ranges of TI values in diagnosis of 
different diseases may be summarized as follows: 

Short TI: High contrast, fat nulled Ti and T 2 weighted sequences for imaging the head 
and neck, abdomen, pelvis, and musculo skeletal systems (e.g., fast STIR or turbo 
STIR). 

Medium TI: Highly Ti weighted sequence with sensitivity to contrast enhancement 
(e.g., MP-RAGE). 

Long TI: CSF-nulled sequence with high sensitivity to acute subarachnoid hemorrhage 
and parenchymal disease, especially disease at the brain-CSF interface, such as multiple 
sclerosis, cortical infarcts, infection, and tumors (e.g., fast FLAIR or turbo FLAIR) 


Description 

Typical 

TI 

values 

Sign of 
longitudinal 
magneti- 
zation 

State of 
longitudinal 
magneti- 
zation 

Usual 
contrast 
weighting 
for target 
tissue 

Tissues 
and fluids 
with Tis 
suitable for 
nulling 

Examples 

Short TI 

80-150 

Negative 

Beginning 

recovery 

Combined 

T 1 &T 2 

Fat, white 
matter, liver 

STIR 

Medium 


Mixed 
negative & 
positive 

Partially 

recovered 

Heavy Ti 

Blood 

(CSF) 

BIR, FIR, 

MP- 

RAGE 

Long TI 


Positive 

Almost fully 
recovered 

Heavy T 2 

CSF 

FLAIR 


Table 1.8: Common clinical variants of the inversion recovery sequence 


Figure 1.6 shows, as given in S. Sinha, M.A. Thomas, A.S. Huda, and R.B. Lufkin 
(2001), the contrast obtained from spin echo images. The images are of a 12 year old 
with HIV and associated CNS lesions. Figure 1.6(a) shows a Ti weighted scan obtained 
using a TE of 14 ms and a TR of 640 ms. Note that Ti contrast dictates that CSF is the 
darkest; white matter is white and gray matter gray. Quite good demarcation is 
observed between the white and gray matter. Figures 1.6(b) and 1.6(c) are obtained 
simultaneously as two images of a double echo sequence. Since a FSE sequence is 
utilized, the relatively short acquisition time permits an increase in the TR to 3800 ms. 























Figure 1.6: Images of a 12 year old with HIV and associated CNS lesions, (a) Ti- 
weighted image with TE/TR: 14 ms/ 640 ms using a conventional spin echo, (b) Proton 
density weighting with a TE/TR; 22ms/3800ms using a Turbo-Spin Echo, (c) A FLAIR 
image with a Turbo-IR, with TE/TR/TI: 105ms/7000ms/2500ms, and (d) A Ta-weighted 
image obtained as the second image in a double echo sequence with the proton density 
image in (b), but with TE: 90ms. Note how the lesion becomes more conspicuous in the 
FLAIR comnared to the T?, and how the CSF is rendered entirely dark. 

The TE for the proton density image is 22 ms and for the T 2 image, 90 ms. Note the 
extremely good quality of the scan, with clear delineation of the different tissues of the 
brain. The contrast in the T 2 weighted scan is opposite compared to the Ti scan, with 



CSF bright, gray matter gray and white matter the darkest. The contrast between tissues 
is lower in the proton density image than in the T 2 image. Finally, in spite of the good 
quality of the T 2 scan, the lesion is almost not visualized. However, when one uses 
FLAIR, the lesion becomes quite distinct. The CSF signal is suppressed but the rest of 
the contrast is similar to T 2 . Very large values of TR and TI are used (7000 ms and 2500 
ms). These would not have been clinically practical but for the use of the Fast Spin Echo 
technique which brings down the scan time to reasonable values with excellent image 
quality. 



1.2.4. Magnetization Transfer Pulse 

Besides the more familiar Ti, T 2 and proton density (p) parameters. Magnetization 
Transfer Contrast (MTC) is one more MR contrast mechanism by means of which 
tissues can be differentiated based on the differences between their intensity values on 
an MR image. Magnetization transfer (MT) utilizes the fact that there is continuous 
magnetization exchange between two proton pools in the brain: the motionally restricted 
pool, which arise from nonaqueous tissue, and the mobile pool from water. MRI can 
only directly detect signal from the mobile pool. If the magnetization from the 
motionally restricted pool is distributed by an MT pulse, then the effect of exchange can 
be seen on the mobile pool as a decrease in signal. This effect is called magnetization 
transfer. This decrease is quantified by means of a magnetization transfer ratio (MTR) 
defined as: 

> f. Me 

MTR= 1 ^ 

I Mo 

where Ms and Mo represent the signals measured with and without off resonance 
radiation respectively. 

Figure 1.7(a) is an image acquired without any MT saturation with TE as 14 ms 
and TR as 1012 ms, and the Figure 1.7(b) is with MT saturation and other sequence 
parameters are same as in Figure 1.7(a). MTR image computed using the above 
equation, shown in Figure 1.7(c), has an excellent contrast between white matter, gray 
matter and CSF. MT effects have started finding increasing clinical usage. Interesting 
use is found in combining MT effect with Gadolinium (Gd) enhancement. MT effects 
do not suppress signal enhancement from Ti shortening caused by Gd infusion in 
extracellular spaces. Consequently, the two effects, parenchymal intensity suppression 


xl00%. 



by MT and lesion enhancement by Gadolinium act synergistically to increase the lesion 
;onspicuity (see J.I. Tanttu, R.E. Sepponen, M.J. Lipton and T. Kuusela (1992)). 



(c) 

Figure 1.7: Magnetization transfer effects. The first image (a) was obtained without MT, the 
second (b) with MT and the third (c) is calculated magnetization transfer ratio, as per equation 
given in text 


1 .3. Literature Review 


In chemistry, computation NMR relaxation times started immediately after the 
discovery of nuclear magnetic resonance phenomena. Many methods and sequences 
were used for computing relaxation times. For example, D.L. De Fontaine, D.K. Ross, 
and B. Temai (1975) described a nonlinear least squares method for the calculation of 
spin-lattice relaxation times. This procedure offered considerable time savings 





compared with other methods, as repetition times only up to the order of Ti are required, 
and, in particular, this procedure obviated the necessity for the determination of S=, the 
signal intensity at thermal equilibrium. 

During the early 1980s, attention was given to parametric mapping since it was 
felt that, in addition to the clinical value of pathology demonstration in parameter- 
weighted images, fundamental information about disease processes and useful 
classification of diseases might be obtained by quantitation of the parameters that yield 
contrast in the weighted images. Attempts were made to obtain this information either 
by pixel-by-pixel calculated parameter maps or by calculations within regions of interest 
in images (M.A. Foster and A. Haase, 1996). Most of the early attempts at Ti and T 2 
mapping used adaptations of the pulse sequences normally employed for in vitro studies, 
e.g., IR for Ti mapping and CPMG, or repeated single spin echo sequences for T 2 . 

A detailed mathematical description of the basic physics and the image 
formation process in MRI could be found in W.S. Hinshaw and A.H. Lent (1983), and 
G. Sebastiani and P. Barone (1991). 

P. A. Bottomley, C. J. Hardy, R. E. Argersinger, and G. Allen-Moore (1987) 
have archived, reviewed, and analyzed the proton NMR relaxation times of pathological 
human and animal tissues in the frequency range 1-100 MHz as a function of tissue of 
origin, NMR frequency, temperature, species, and in vivo versus in vitro status. Tj data 
from specific disease states of the bone, brain, breast, kidney, liver, muscle, pancreas, 
and spleen was characterized by simple dispersions of the form Ti = Av® in the range 1- 
100 MHz with A and B empirically determined pathology-dependent constants. 
Pathological tissue T 2 values are essentially independent of NMR frequency. Raw 
relaxation data, best-fit Tj parameters A and B, and the mean T 2 values, were tabulated 
along with standard deviations and sample size to establish the normal range of 
pathological tissue relaxation times applicable to NMR imaging or in vitro NMR 
examination. They have observed from statistical analysis of relaxation data that most 
tumor and edematous tissue Ti values and some breast, liver, and muscle tumor T 2 
values were significantly elevated (p > 0.95) relative to normal, but did not differ 
significantly from other tumors and pathologies. Statistically significant abnormalities in 
the Ti values of some brain, breast, and lung tumors, and most pathological tissue T 2 
values could not, however, be demonstrated in the presence of large statistical errors. 
Both Ti and T 2 in uninvolved tissue from tumor-bearing animals or organs did not 
demonstrate statistically significant differences from normal when considered as a 
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group, suggesting no appreciable systemic effects associated with the presence of 
tumors compared to the statistical uncertainty. Statistical prediction analysis for both Ti 
and T 2 indicated that of all the tissues studied, only liver hepatoma could be reliably 
distinguished from normal liver based on a single T 1 measurement (p > 0.95) given the 
scatter in the published data. Indeed, data scatter, not easily attributable to temperature, 
species, in vivo versus in vitro status, the inclusion of implanted or chemical induced 
tumors, or the possible existence of multiple component relaxation, was recognized as 
the major factor inhibiting the diagnostic utility of quantitative NMR relaxation 
measurements. 

L.R. Schad, G. Brix, I. Zuna, W. Harle, W.J. Lorenz, and W. Semmler (1989) 
have performed the in vivo measurements of proton relaxation processes in human brain 
tumors by MR imaging using a whole-body superconductive MR scanner, operating at 
1.5 T. The Ti and T 2 relaxation time measurements were based on a combined Carr- 
Purcell/Carr-Purcell-Meiboom-Gill sequence with two interleaved repetition times and 
32 echoes. First, comparative measurements in the imager and with spectrometer of 
relaxation times were performed on a phantoms containing fluids of different Ti and T 2 
to evaluate accuracy. A maximum deviation of ~10% was found. Multislicing with a 
gap width of one slice thickness influenced the accuracy of T i relaxation measurement. 
A gap width of at least two times the slice thickness was necessary for reliable 
determination of Ti. No influence on T 2 values was observed by multislicing. Second, in 
human head imaging the multiexponential behavior of the T 2 decay curves was analyzed 
in each pixel, where the mean square deviation were used as a creterian to discriminate 
between mono and biexponential behavior. Mean values of monoexponential Ti and 
multiexponential T 2 relaxation data for white matter, gray matter, CSF, edema, and 
tumor were sampled in 12 patients with brain tumors. These authors observed that T 2 
showed mono exponential behavior in white matter and gray matter, whereas CSF, 
edema, and tumor showed distinct biexponentiality. The biexponential analysis 
generally yields "fast" and "slow" components with T 2 f = 80 ± 17 ms and T 2 S = 2,030 ± 
210 ms for CSF (partial volume effect), T 2 f = 104 ± 25 ms and T 2 s = 677 ± 152 ms for 
edematous tissues, T 2 f = 97 ± 19 ms and T 2 s = 756 ± 99 ms for tumor tissues, 
respectively. For the discrimination of edematous and tumor tissues a retrospective 
overall accuracy of 94% was found. 

H. Iwaoka, T. Hirata and H. Mastuura (1987) have tested several scan sequences. 
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of Ti, T 2 , and proton density values encountered in representative human tissues. They 
found that, for given scan time, the combination of inversion recovery 3 spin echo 
(IR3SE) and saturation recovery 4 spin echo (SR4SE) sequences gave more accurate 
computed images than other comparable methods tested. 

J. Gong and J.P. Hornak (1992) have addressed an improvement in the post 
processing time required to generate Tj images. They have described a nonlinear least- 
squares algorithm for rapidly generating spin-lattice relaxation time images from 
variable repetition time magnetic resonance images. The algorithm generates a 256x256 
pixel Ti image from nine variable repetition time images in approximately 60 sec. on a 
VAX-6510 computer. 

Xing Li and J.P. Hornak (1994) have generated T 2 images are from a series of 
variable echo-time magnetic resonance images by either a linear or a nonlinear least- 
squares technique. The linear technique is contrasted with a rapid nonlinear technique 
developed by the authors. Compared with linear techniques, the nonlinear technique 
produced calculated T 2 images with less noise and less error in those tissues with low T 2 
values, such as adipose and muscle. 

Calculated Ti images require that magnetic resonance signal be detected as 
several inversion or repetition times. Multiple spin-echo acquisitions provide several 
measurements of the magnetization at each TR, the signal size diminishing according to 
T 2 decay. S.J. Riederer, S.A. Bobman, J.N. Lee, F. Farzaneh, and H.Z. Wang (1986) 
have reviewed one method (case 1) for estimating Ti from single echoes and presented 
four new methods (case 2-5) in which multiple acquired echoes were used. For case 2 a 
fit was performed using the first echo at each TR, repeated using second echo, etc., and 
the final Ti estimate was the simple average of the individual fits at each echo time 
(TE). For case 3 the optimum weighted average was performed. For cases 4 and 5 
synthetic SE images were generated at each TR prior to the Ti fit, case 4 using a 
synthetic TE of zero, and case 5 using a TE providing maximum signal-to-noise ratio in 
the synthetic image. The relative precision in Ti provided by each method was 
calculated. It was proven that case 3 and 5 are optimum and equivalent and could 
theoretically reduce the noise in Ti images by as much as 40% over case 1 with no 
increase in scanning time. Approximations were proposed that enable the optimum 
methods to be implemented in a practical fashion. 

Techniques used in NMR imaging, not routinely employed in conventional 
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For the saturation recovery approach of Ti measurement, these include selective 
excitation to define the image plane and magnetization refocusing for NMR signal 
acquisition. B.R. Rosen, L.I. Pykett, and T.J. Brady (1984) showed by computer 
modeling of the Bloch equations that errors well over 50% could be made in image 
derived Ti measurements if the conventional SR relations between signal intensity and 
the 90-90° inter pulse delay, TR, is used. They have derived corrected expressions for 
the actual pulse and gradient sequence used by their imaging system, and phantom data 
acquired in imaging experiments have verified the validity of these equations. This 
allowed for the correction of Ti data to compensate for systematic bias introduced 
during imaging procedures and suggested a means whereby errors introduced by radio 
frequency inhomogeneities across the imaging volume could be reduced. 

I. Harvey, P.S. Tofts, J.K. Moris, D.A.G. Wicks, and M.A. Ron (1991) have 
studied the major factors contributing to Tj variance at 0.5 T in white matter in healthy 
people. Anatomical location of the white matter sampled and differences between 
individuals contributed 74% of the total variance in serial measurements of the same 
objects. There was also significant change over time within an individual subject that 
could not be attributed to machine drift. This information permitted estimates to be 
made concerning adequate sample size in studies that examine for pathological white 
matter Ti change. 

Some of the major problems encountered when measuring relaxation times in vivo 
are (M.A. Foster and A. Haase, 1996): 

1 . magnetic field inhomogeneity effects, including those due to the field gradients, 

2. slice profile effects, related to pulse sequence timing etc., 

3. problems relating to rf pulses and signal phase, 

4. errors when multislice collection procedures are used, 

5. misregistration problems when Ti and T 2 are calculated from several imaging data 
sets, 

6. partial volume effects due to tissue inhomogeneity, 

7. assumptions made in the Ti and T 2 calculation, e.g., regarding relaxation 
monoexponentiality, 

8. a wide variety of effects arising from the living tissue/body itself, such as flow, 
diffusion, magnetic susceptibility, etc. which can all either affect the relaxation time 
value directly or alter the apparent place from which the information is coming. 
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S.A. Bobman, SJ. Riederer, J.N. Lee, S.A. Suddarth, H.Z. Wang, J.R. MacFall 
(1985) have studied the comparison of synthesized MR images with acquired images. 
They have found that the synthetic images subjectively compared favorably in both 
accuracy and precision with acquired images when formed for the same values of echo 
and repetition times and for interpolated and extrapolated values of both TE and TR. 
Percent error was consistently less 5% for brain parenchyma, and synthetic signals were 
accurate to within four times the noise level at acquisition. The apparent signal-to-noise 
ratio of synthetic images was comparable, superior, or inferior to similar acquired 
images, depending on the values of TE and TR. S.A. Bobman, S.J. Riederer, J.N. Lee, 
T. Tasciyan, F. Farzaneh, H.Z. Wang (1986) have done the pulse sequence extrapolation 
with MR image synthesis. They have used the multiple spin-echo data acquired in a time 
equivalent to a TR of 2,000 ms and generated inversion-recovery images for arbitrarily 
chosen inversion times. Synthetic IR signals of the brain parenchyma consistently 
matched directly acquired signals to within 6%. 

Any chemical exchange or magnetization transfer between the water and 
macromolecular protons will result in a decease in the water proton signal or saturation 
transfer as originally described by S. Forsen and R.A. Hoffman (1966) for simple 
chemical and dipolar interactions. The demonstration of normally unobservable proton 
pools in exchange with the observable pools is an important finding which increases the 
potential image contrast parameters available in MR imaging. S.D. Wolff and R.S. 
Balaban (1989) first produced in vivo images with magnetization transfer weighted 
contrast. They have studied the exchange between ^H magnetization in "free" water 
('Hf) and that in a pool with restricted motion ('Hr) observed in tissues in vivo using 
NMR saturation transfer methods. Exchange between these two pools was demonstrated 
by a decrease in the steady-state magnetization and relaxation times of 'Hf with 
radiofrequency irradiation of 'Hr. Proton NMR imaging demonstrated that this exchange 
was tissue specific and generated a novel form of NMR image contrast. 

B.S. Hu, S.M. Conolly, G.A. Wright, D.G. Nishimura, and A. Macovski (1992) 
have examined the theoretical behavior of short Ta species in the presence of RF 
irradiation with regard to pulse duration, amplitude, and total integrated RF power, as 
well as off-resonance effects. J.V. Hajnal, CJ. Baudouin, A. Oatridge, I.R. Young, and 
G.M. Bydder (1992) have discussed the design and implementation of magnetization 
transfer pulse sequences for clinical use. 
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A good amount of literature on the applicability of magnetization transfer effect 
and the quantitative magnetization transfer ratios for medical diagnosis of different 
diseases is available. R.C. Mehta, G.B. Pike, D. Enzmann (1996) have studied the 
measure magnetization transfer in multiple sclerosis demyelinating plaques, white 
matter ischemic lesions, and edema. J.R. Petrella, R.I. Grossman, J.C. McGowan, G. 
Campbell, J.A. Cohen (1996) have analyzed the relationship between MR enhancement 
pattern and magnetization transfer effect in multiple sclerosis lesions. There are a 
number of studies e.g., A. Gass, G.J. Barker, D. Kidd, et al. (1994), J.F. Hiehle, R.I. 
Grossman, K.N. Ramer, F. Gonzalez-Scarano, J.A. Cohen (1995), L.A. Loevner, R.I. 
Grossman, J.C. McGowan, K.N. Ramer, J.A. Cohen (1995), M. Fillippi, M.D. Campi, 
V. Dousset, et al. (1995) carried out on the magnetization transfer effect on normal 
appearing white matter in multiple sclerosis. 

Measurement of magnetization transfer in different stages of 
neurocysticercosis was carried out by M.K. Kathuria, R.K. Gupta, R.Roy, V. Gour, N. 
Hussain, and S. Pradhan (1998). R.K. Gupta, M.K. Kathuria and S. Pradhan (1999) have 
found that an MT effect increases the detectability of CNS tuberculosis and helps in 
better assessment of the disease load. 

1 .4. Contents of the Thesis 

A study of the computational methods for tissue parameters estimation and the 
effect of perturbations in signals on parameter estimation is discussed in chapter 2. 

In section 2.1, most commonly used two point methods for estimating Ti and T 2 
are presented. Errors due to the choices of TR's for Ti estimation and TE's for T 2 
estimation in the presence of noise for some brain tissues are presented. 

In section 2.2, some iterative methods for tissue parameter estimation are 
described. Subsection 2.2.1 describes a least squares method given by J. Gong and J.P. 
Hornak (1992) for computation of Ti. This method requires the data acquired at 
different TR's keeping TE fixed. A look-up table approach is suggested for faster 
computation of Ti values compared to the controlled iterative process suggested by the 
authors. The accuracy in the computed Ti values is improved by Newton's method by 
using the look-up table values as initial approximations. In subsection 2.2.2, a nonlinear 
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constant. For this method also a look-up table approach is suggested to make this 
method faster than the controlled iterative process proposed by the authors. The 
accuracy in the computed T 2 values is improved by Newton's method by using the look- 
up table values as initial approximations. Subsection 2.2.3 describes a least squares 
method given by H. Iwaoka, T. Hirata and H. Matsuura (1987) for the computation of 
Ti, T 2 , and PD values simultaneously. Given three or more independent images 
corresponding to a cross-section, this method computes Tj, T 2 , and PD values 
iteratively. The convergence of this method depends on the initial approximations for 
Ti, T 2 and PD values. In practice, these initial approximations are not known before 
hand. Any look-up table kind of approaches would take a lot of time in searching for 
these initial values due to the high degree of freedom available in this method and the 
memory space required. In subsection 2.2.4, we have presented a new method which 
reduces the system of equations to one involving a single function of single variable 
case. Given independent data set of three images, it is shown by proper re-arrangement 
of the signal expressions that one could reduce them in to a function of Ti. Then Ti 
values are obtained from a generated look-up table. Using these as initial 
approximations, the accuracy in the computed values could be improved iteratively by 
Newton's method. By substituting these Ti values into any two image intensities, T 2 
could be computed. This T 2 computation does not require any specific data, unlike in the 
method described in subsection 2.1.2. Now, PD values are obtained by directly 
substituting the computed Ti and T 2 values in to any one of the image intensity 
equations. If more than three images are available, by considering suitable data triples, a 
least-squares method is suggested to compute more accurate Ti values and consequently 
T 2 and PD values. Subsection 2.2.5 presents a least squares method for computing Ti, T 2 
and PD values simultaneously. Initial approximations for this iterative method are 
obtained from the computed values using the method described in subsection 2.2.4. 

In all the methods described above, the signal amplification proportionality 
constant is assumed to be constant. But in practice, this constant may vary over the 
regions or images due to gradient and signal receiver coil imperfections. In subsection 
2.2.6 a method is presented for the estimation of tissue parameter computables like Ti, 
Tlx, T 2 , T 2 x, P, and Px in the presence of amplification factors where the suffix x denotes 
the variation in x-direction. Computation of the field susceptibility maps is presented. 
Simulation results are given to demonstrate this method. 
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A numerical study of the sensitivity of the tissue parameters computation due to 
perturbations in the input images is presented in section 2.3. In subsection 2.3.1, errors 
in the computed parameters due to bit reduction in the images are presented. Subsection 

2.3.1. deals with the errors in parameters estimation due to blurring in the input images. 

Chapter 3 mainly deals with the generation of high contrast images. In section 

3.1, an entropy-like maximization procedure which computes the weights is presented. 
Using these weights, a linear combination of the Ti, T 2 , and PD weighted acquired 
images is produced. The resultant image is similar to PD weighted image, but has shown 
a good contrast between perilesional gliosis region and brain parenchyma than in any 
individual image. Given a set of linearly independent images, Gram-Schmidt 
orthogonalization process is discussed in section 3.2 to generate a set of orthonormal 
images. Results of orthonormalizing Ti weighted image with respect to T 2 weighted 
image and T 2 weighted image with respect to Ti weighted image are presented. In 
section 3.3, generation of contrast images through functional maps is discussed. By 
substituting the computed tissue parameter values into the functions of Ti and T 2 (e.g., 
for saturation recovery spin echo, /(Ti) = l-2exp(-(TR-TE/2)/Ti)+exp(-TR/Ti), and, 
g(T 2 ) = exp(-TE/T 2 )) at desired TE and TR values, an excellent contrast in the images is 
produced. In section 3.4, a least squares method is presented to produce tissue parameter 
maps contrast through a weighted combination of the images. The weights computed at 
one slice did not vary much from the weights associated with adjacent slices. So, 
weights computed at one slice could be used to produce tissue parameter maps contrast 
like images without computing these maps. In section 3.5, generation of good contrast in 
perilesional gliosis region on synthesized IR images is demonstrated. Synthesis of image 
values at divergent point, where tissue parameter values could not be computed due to 
partial volume effect, is discussed in section 3.6. It is observed that the tissue parameter 
values in different segments (regions) of brain are normally distributed. Utility of tissue 
parameter maps for segmentation using a variant of MLE algorithm (R.K.S. Rathore, S. 
Datta, R.K. Gupta, S.B. Rao and R. Verma (2001)) is shown in section 3.7. 

Chapter 4 is dedicated to the quantitative analysis of the regions of three 
pathologies: cystic cavity, perilesional gliosis and edema. In section 4.1, the statistical 
tools used for the analysis are discussed. The means and variances in pathological 
regions and contra-lateral normal regions on MT SE Ti, Ti, T 2 and PD weighted images 
as well as on MTR, Ti, T 2 and PD computed maps have been analyzed. In all these 
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the significance of the changes in the means and variances. The linear correlation 
coefficients between different weighted images and between computed maps in normal 
and pathological regions are presented. Section 4.2 is devoted to a quantitative analysis 
of cystic cavity region. On MT SE Ti, Tj and PD weighted images the means in the 
cystic cavity regions are less compared with contra-lateral normal regions. But T 2 
weighted images have shown high values in the lesion region in comparison with 
normal region. Ti, T 2 and PD maps have shown high mean values, while MTR map has 
less values in the lesion regions compared with normal regions. In section 4.3, 
perilesional gliosis around the cystic cavity has been studied. Analysis suggests the 
usefulness of magnetization transfer for gliosis. On Ti maps, the mean values in the 
lesion are significantly high compared with normal region even though it is not observed 
on T 1 weighted images. Analysis of edema regions is discussed in section 4.4. 



Chapter 2 

Numerical Methods for Tissue Parameter 
Quantitation and Sensitivity Analysis 


Knowledge of Ti and T 2 is not simply an academic curiosity. Many diseases are 
apparent with certain pulse sequences but not visible with others (e.g., multiple sclerosis 
plaques). The fact that these lesions can be seen at all shows that they have different p, 
Ti, and/or T 2 values from normal tissues, but the fact that they are not seen on some 
sequences means that an optimization scheme is needed to assure that the proper 
sequence is used without spending excessive time to find it. In one such scheme, M.R. 
Mitchell, T.E. Conturo, T.J. Gruber, and J.P. Jones (1984) are able to optimize the 
contrast based upon a knowledge of Ti and T 2 . Further, there is a possibility that p, Ti, 
and T 2 may aid in the differential diagnosis of disease. Without knowledge of T] and T 2 , 
MR imaging will remain a trail and error technique with good sensitivity but low 
specificity. In other words, some of its potential may be wasted (J.P. Jones (1988b)). 

A study of the computational methods for tissue parameters estimation and the 
effect of perturbations in signals on parameters is discussed in this chapter. 

In section 2.1, most commonly used two point methods for estimating Ti and T 2 
are presented. Errors due to the choices of TR's for Ti estimation and TE’s for T 2 
estimation in the presence of noise for some brain tissues are presented. 

In section 2.2, some iterative methods for tissue parameter estimation are 
described. Subsection 2.2.1 describes a least squares method given by J. Gong and J.P. 
Hornak (1992) for computation of Ti. In subsection 2.2.2, a nonlinear least squares 
method for computation of T 2 values given by Xing Li and J.P. Hornak (1994) is 
described. Subsection 2.2.3 describes a least squares method given by H. Iwaoka, T. 
Hirata and H. Matsuura (1987) for the computation of Ti, T 2 and PD values 
simultaneously. In subsection 2.2.4, we have presented a new method which reduces the 
system of equations to one involving a single function of single variable case. A new 
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where suffix x denotes the variation in x-direction, in the presence of amplification 
factors is presented in subsection 2.2.6. 

A numerical study of the sensitivity of the tissue parameters computation due to 
perturbations in the input images is presented in section 2.3. In subsection 2.3.1, errors 
in the computed parameters due to bit reduction in the images are presented. Subsection 
2.3.1 deals with the errors in parameters estimation due to blurring in the input images. 


2.1. Two Point Calculations 

Both the magnetic resonance signal intensity and contrast between tissues 
depend upon the local values of spin density (p), and the relationships between Ti, T 2 , 
and the interpulse time intervals. Neglecting flow, signal expression can be written as 

S = kpf(T,)g{T^), 

where S is the signal intensity, ^ is a proportionality constant, and/& g are functions of 
the type of pulse sequence (J.P. Jones (1988)). In the case of spin echo sequence,/ and 
g are given by: 

/ = 1 + exp(-ri? / r, ) - 2 exp(-(r/? -TE/2)/Ti), 
g =Qxp{-TE/T2), 

where TR and TE are the repetition and echo times respectively. 

Calculating Tp Keeping TE fixed while collecting images with different TR values, 
the ratio of signal intensities is 

5, (1 + expC-TRi / Ti ) - 2 exp(-(ri?i -TEH)IT ^ )) 

^2 ” (1 + expC-TRj / Ti ) - 2 expC-CTRj - TE / 2) / / )) ’ 

This can be solved for Ti. In practice, TR»TE/2, so this ratio can be written as 

Si_ ^ (l-expl-TRi/Ti)) ^ 

^2 (l-exp(-77?2/T,))’ 

and the value of Ti can be deduced by an iterative scheme. Initial guesses could be 
obtained from a lookup table. In addition, if TRi=2TR2, from the above equation it can 
be written as 


Ti =-ri?2/ln((5,/52)-l). 


(2.1.2) 



Figure 2.1(a) and Figure 2.1(b) are the two single spin echo images with same 
TE = 14 ms and with different TR's, 1000 ms and 2000 ms respectively. Figure 2.1(c) is 
the corresponding Ti map computed using the equation (2.1.2). 



(a) TE/TR = 20/2200 ms 


(b)TE/TR = 20/1 100 ms 



(c)Ti map 

Figure 2.1: (a) and (b) are two spin echo images with TE/TR = 20/2200, 1100 ms. Image (c) is 
the computed Ti map 


Percentage errors are computed in Ti estimations of some common brain tissues 
using equation (2.1.2) with different TR's at different noise levels. First, at each brain 
tissue Ti value, for given two TR values the true signal ratio S 1 /S 2 as in equation (2.1.1) 
is computed. Different percentages of these ratios are added as noise and resulted signal 
ratios are used in equation (2.1.2) to compute Ti values (J.P. Jones (1988b)). Tables 2.1- 
2.4 show the percentage errors in computed Tivalues of different brain tissues in the 





presence of different noise levels in the signal ratios using the equation (2.1.2) with 
(TRi, TR 2 ) = (1200, 600), (1500, 750), (2000, 1000) and (3000, 1500) ms respectively. 
Table 2.1 shows at (TRi, TR 2 ) = (1200, 600) ms the estimation of Ti's of brain tissues 
with high Tivalues are more sensitive to the presence of nose compared to the tissues 
with low Ti's. For example, at 1% noise methemoglobin and CSF have shown errors 
3,64 and 10.92 respectively, and at 10% noise while methomoglobin has error 41.76, 
CSF increased up to 878.46. From Tables 2. 1-2.4, it is observed that at high range of 
(TRi, TR 2 ) values, errors in computation of Ti values of tissue with high Tivalues like 
CSF, upper cystic fluid, edema, gray matter and lower cystic fluid are less compared to 
the errors in estimation at lower (TRi, TR 2 ) values. But tissues with low Ti values like 
methemoglobin and white matter have shown increased errors at higher range of (TRi, 
TR 2 ) values. 


Tissue Type 

Ti 

At 1 % Noise 

At 5% Noise 

At 10% Noise 

Ti 

%Error 

T, 

%Error 

Ti 

%Error 

White matter 

510 

528.67 

3.66 

609.71 

19.55 

729.22 

42.98 

Gray matter 

760 

791.60 

4.16 

936.10 

23.17 

1172.62 

54.29 

Edema 

900 

941.01 

4.56 

1133.83 

25.98 

1469.40 

63.27 

CSF 

2650 

2939.38 



89.28 

25929.21 

878.46 

Upper Cystic Fluid 

1080 

1135.30 

5.12 

1405.04 

30.10 

1915.92 

77.40 

Lower Cystic Fluid 

720 

749.20 

4.05 

881.61 

22.45 

1094.73 

52.05 

Methomoglobin 

460 

476.73 

3.64 

548.51 

19.24 

652.11 

41.76 

Table 2.1: % Errors in Ti computation of 
= 1200 ms and TR 2 = 600 ms at different 1 

'different brain tissues using equation (2.1.2) with TRi 
loise levels 

Tissue Type 

■ 

At 1% Noise 

At 5% Noise 

At 10% Noise 

Ti 

%Error 

T, 

%Error 

Ti 

%Error 

White matter 

510 

528.74 

3.67 

608.04 


719.81 

41.14 

Gray matter 

760 

788.91 

3.80 

917.06 

20.67 

1113.93 

46.57 

Edema 

900 

936.50 

4.05 



1368.41 

52.05 

CSF 


2884.46 

8.85 

4336.55 

63.64 


283.49 

Upper Cystic Fluid 

1080 

1128.06 

4.45 



1736.59 

60.79 

Lower Cystic Fluid 

720 

746.98 

3.75 

865.76 

20.24 

1045.84 

45.25 

Methomoglobin 

460 

477.35 

3.77 

549.85 

19.53 

650.01 

41.31 


Table 2.2: % Errors in Ti computation of different brain tissues using equation (2.1.2) with TRi 
= 1500 ms and TR 2 = 750 ms at different noise'levels 






























































































Tissue Type 

wnm 

At 1 % Noise 

At 5% Noise 

At 10% Noise 


m 

T, 

%Error 

T. 


T, 

%Error 

White matter 

510 

531.11 

4.14 





Gray matter 

760 

787.65 

3.64 





Edema 

900 

933.25 

3.69 

1078.49 

19.83 



CSF 




3825.19 

44.35 

6347.34 

139.52 

Upper Cystic Fluid 

1080 

1121.97 


1309.54 

21.25 

1602.49 

48.38 

Lower Cystic Fluid 

720 

746.27 

3.65 

858.13 

19.18 

1017.57 

41.33 

Methomoglobin 

460 

480.66 

4.49 

563.26 

22.45 

670.61 

45.78 

Table 2.3: % Errors in Tj computation of different 1 
= 2000 ms and TRo = 1000 ms at different noise leve 

jrain tissues using equation (2.1.2) with TRj 
;ls 

Tissue Type 

D 

At 1 % Noise 

At 5% Noise 

At 10% Noise 



T, 


T, 

%Error 

Ti 

%Error 

White matter 

. 

543.60 

6.59 

666.79 

30.74 



Gray matter 

760 

791.60 






Edema 

900 

934.22 

3.80 

1076.78 

19.64 

1272.89 

41.43 

CSF 








Upper Cystic Fluid 

1080 






H 

Lower Cystic Fluid 

720 

751.18 

4.33 

876.84 

21.78 



Methomoglobin 

460 

496.48 

7.93 

623.72 

35.59 

769.00 

67.17 


Table 2.4: % Errors in Tj computation of different brain tissues using equation (2.1.2) with TRi 
= 3000 ms and TR 2 = 1500 ms at different noise levels 


As TRi and TR 2 are related (TRi=2TR2), Graphs 2.1 (a )-2.1(g) show the % 
errors in Ti computation against different TR 2 ( = 600, 750, 1000,1500 ms) values in 
white matter, gray matter, edema, CSF, upper cystic fluid, lower cystic fluid and 
methemoglobin respectively at 1%, 5% and 10% noise levels. Curves are plotted with 
spline interpolation. Graphs 2.1(a) and 2.1(f) show that gray matter and LCF have less 
errors at TR 2 =1000 ms compared to the other TR 2 values. Graph 2.1(c)-2.1(e) and 
2. 1 (g) show the decrease in errors with increase in TR 2 values. 















































































% Errors in computation of Edema % Errors in computation of GM %Erras in computation of WM 



Graph 2.1 (a): % Errors in Ti 
different noise levels 


computation of white matter with different TR 2 values at 



TRj Values 


Graph 2.1 (b): % Errors in Ti computation of gray matter with different TR 2 values at 
different noise levels 



Graph 2.1 (c): % Errors in Ti computation of edema with different TR 2 values at 
different noise levels 


%Errors in computation of LCF %Errors in computation of UCF 



TRj Values 

Graph 2.1 (d); % Errors in Ti computation of CSF with different TR 2 values at 
different noise levels 



TRj Values 


Graph 2.1 (e): % Errors in Ti computation of upper cystic fluid with different TR 2 
values at different noise levels 



Graph 2.1 (f): % Errors in Ti computation of lower cystic fluid with different TR 2 values 
at different noise levels 




Graph 2.1 (g): % Errors in Tj computation of methemoglobin with different TR 2 values 
at different noise levels 


Calculating T2: Using the approximation TR»TE/2, the value of T2 can be deduced by 
holding TR constant while varying TE. The ratio is then 

5i _ expC-TEi /r^) 

^3 exp(-T£'3 /Tj) ’ 

from which T2 directly follows as; 

T2=iTEj-TE^)/\n(S^/S^). 

Figure 2.2(a) and Figure 2.2(b) are, respectively, PD and T2 weighted spin echo 
images obtained from a double echo sequence with TR as 2200 ms and TEi as 12 ms 
and TE2 as 80 ms. Figure 2.2(c) shows the T2 map computed as given in the equation 
(2.1.4). 

Percentage errors are computed in T2 estimations of some common brain tissues 
using equation (2.1.4) with different TE's at different noise levels. First, at each brain 
tissue T2 value, for given two TE values the true signal ratio S2/S3 as in equation (2.1.3) 
is computed. Different percentages of these ratios are added as noise and resulted signal 
ratios are used in equation (2.1.4) to compute T2 values (J.P. Jones (1988b)). Tables 2.5- 
2.8 show the percentage errors in computed T2 values of different brain tissues in the 
presence of different noise levels in the signal ratios using the equation (2.1.2) with 
(TEi, TE2) = (20, 60), (20, 80), (20, 140) and (20, 200) ms respectively. From these 
tables it is observed that the errors in T2 computation for all tissues at different noise 
levels decrease as the difference between TEi and TE2 increases. In acquired MR 


(2.1.3) 

(2.1.4) 



images the SNR will be low at higher TE values and consequently T 2 values computed 
with low TEi and high TE 2 correspond to the case of the presence of high noise level. 
For example, from Tables 2.7 and 2.8, the T 2 values computed for CSF with (TEi, TE 2 ) 
= (20, 200) at 10% noise level has an error 12.91 while with (TEi, TE 2 ) = (20, 140) at 
5% noise level has an error 10.22 only. 



(c) T 2 map 


Figure 2.2: (a) and (b) are two spin echo images with TR/TE = 2200/20,80 ms. Image (c) is 
the computed T 2 map 





Tissue Type 

■ 

At 1% Noise 

At 5% Noise 

At 10% Noise 

T 2 


Ta 

%Error 

T 2 


White matter 

67 

65.90 

1.64 

61.94 

7.55 

57.78 

13.77 

Gray matter 

77 

75.55 

1.88 

70.39 

8.59 

65.06 

15.50 

Edema 

126 

122.17 

3.04 

109.21 

13.32 

96.91 

23.09 

CSF 



6.51 

208.72 

25.46 

167.95 

40.02 

Upper Cystic Fluid 






142.17 

33.87 

Lower Cystic Fluid 

43 



40.86 

4.98 

39.00 

9.29 

Methomoglobin 

106 



93.86 

11.45 

84.63 

20.16 


= 20 ms and TE 2 = 60 ms at different noise levels 


Tissue Type 

■ 

At 1% Noise 

At 5 % Noise 

At 10% Noise 

T 2 

%Error 

T 2 

%Error 

T 2 


White matter 

67 

66.26 


63.54 

5.12 

60.55 

9.62 

Gray matter 

77 

76.03 


72.46 

5.89 

68.61 


Edema 

126 

123.42 

2.05 

1 14.29 

9.29 

104.99 

16.68 

CSF 




228.07 

18.55 

193.80 

30.78 

Upper Cystic Fluid 







25.46 

Lower Cystic Fluid 

43 

42.69 

0.71 

41.55 

3.38 


6.39 

Methomoglobin 



1.73 

97.59 

7.94 


14.41 


Table 2.6: % Errors in T 2 computation of different brain tissues using equation (2. 
= 20 ms and TE 2 = 80 ms at different noise levels 


.4) with TEi 


Tissue Type 

Tz 

At 1 % Noise 

At 5% Noise 



At 10% Noise 



Tz 






White matter 

67 

66.63 

0.55 

65.22 




Gray matter 




74.66 

3.04 

72.56 

5.76 

Edema 

126 

124.70 

1.03 

119.86 

4.87 



CSF 

280 



251.38 

10.22 

229.06 

18.19 

Upper Cystic Fluid 



1.75 



183.64 

14.59 

Lower Cystic Fluid 

43 

42.85 

0.35 

42.26 

1.72 

41.58 

3.30 

Methomoglobin 






97.77 

7.76 


= 20 ms and TE 2 = 140 ms at different noise levels 




























































































Tissue Type 

■ 


At 5% Noise 

At 10% Noise 

T 2 

%Error 

T 2 

%Error 

T 2 


White matter 

67 

66.75 

0.37 

65.80 

1.78 

64.70 

3.43 

Gray matter 

77 

76.67 

0.42 

75.43 

2.04 


3.92 

Edema 

126 

125.13 

0.69 

121.84 

3.30 

118.12 

6.25 

CSF 

280 

275.73 



7.05 

243.85 

12.91 



212.47 

1.17 

203.16 

5.51 



Lower Cystic Fluid 

43 

42.90 

0.24 

42.50 

1.15 



Methomoglobin 



0.58 

103.04 

2.79 




Table 2.8: % Errors in T 2 computation of different brain tissues using equation (2.1.4) with TEi 
= 20 ms and TE 2 = 200 ms at different noise levels 


As TE| is 20 ms in all the (TEi, TE 2 ) pairs. Graphs 2.2(a )-2.2(g) show the % 
errors in T 2 computation against different TE 2 ( = 60, 80, 140, 200 ms) values in white 
matter, gray matter, edema, CSF, upper cystic fluid, lower cystic fluid and 
methemoglobin respectively at 1%, 5% and 10% noise levels. Curves are plotted with 
spline interpolation. All these graphs show that the errors in T 2 computation decreases 
with the increase in TE 2 values. 



Graph 2.2 (a): % Errors in T 2 computation of white matter with different TE 2 values at different 
noise levels 









































Graph 2.2 (b): % Errors in To computation of g 
noise levels 



Graph 2.2 (c): % Errors in T 2 computation ol 
levels 





Graph 2.2 (e): % Errors in T 2 computation of upper cystic fluid with different TE 2 values at 
different noise levels 



Graph 2.2 (f): % Errors in T 2 computation of lower cystic fluid with different TE 2 values at 
different noise levels 



Graph 2.2 (g): % Errors in T 2 computation of methemoglobin with different TE 2 values at 
different noise levels 


2.2. Iterative Methods 


Usually, nonlinear equations are solved numerically by iterative methods. Tissue 
parameter estimation in MRI is a nonlinear problem in nature. In this section we discuss 
some iterative methods for computation these parameters. 


2.2.1. A Least Squares Method for Ti Calculation 


In this subsection we describe a least squares method proposed by J. Gong and J.P. 
Hornak (1992) for the computation of Ti values. Under ideal conditions, that is, the 
absence of spin exchange, primarily one component filling an image voxel, and perfect 
radio-frequency pulses, the magnetic resonance signal from a saturation recovery type 
imaging experiment may be approximated as 

S{TRj ) = k(l- expi-TRj /T ^ )) (2.2.1) 

where TR is the repetition time of the experiment and ic is a proportionality constant 
dependent on such factors as the spin density and amplifier gains. In the case of 
commonly used spin-echo sequence, k is also a function of the spin-spin relaxation 
time, TE. A Ti value is typically calculated by fitting equation (2.2.1) to imaging data 
obtained at several values of TR. The imaging data is the actual MR signal Si at TRi. A 
least squares technique seeks to minimize (j)(Ti, k) for the n data points, where 

m^k) = '^[S(TR,)-Sf. 


Therefore 


S<l>(T,,k)=^Sk+^gr, = o. 


Since neither 6k nor 6Ti may be zero then 

\j aiiu ^ 

dk dT^ 


0 and ^ = 0. 


From the above equation one has 

^[y(:(l - exp(-TRi ITJ}- S, ](l - exp(-TR,- /Tj )) = 0 
1=1 


and 


£ [ife(l - exp(-rR, IT,))- 5, FR, exp(-rR, / Tj ) = 0 . 



X5,.(l-exp(-r/?, IT,) exp(-r/?,. IT,) 

JM - tl 

^ (1 - txp{-TR, IT,)^ '^TR, exp(-77?,. / 7, )(1 - exp(-Ti?,. / T, )) 

i=l ^'=1 

Now the above equation is a function of Ti only. After reordering the summations one 
obtains a function Z(Ti) given as 

^(^i) = IT,){l-txpi-TR, IT,) 

i=i j=\ 


X [S; (1 - ex^i-TRj /T,))-Sj{l- exp{-TR, / T , ))] } . 


The authors have used a controlled iterative method to computed Ti values. 

This least squares method requires the data acquired at different TR values but 
with same TE values. It is also to be noted that the signal expression used in this method 
is an approximated expression i.e. with the assumption TR » TE. 

Figure 2.3(a)-(d) shows the spin echo images of a volunteer acquired with 
TR/TE = 850, 1100, 1650, 2200/20 ms. Figure 2.3(e) is the Tj map computed from 
Figure 2.3(a)-(d) images using the above method. White dots outside the brain Ti map 
are due the presence of noise more than the threshold value, here it is taken as 60, in the 
input images. Computing T i map with this data using the method described above has 
taken 16.04 sec on a Pentium II P.C. with 350 MHz speed and 96 MB RAM. 

This method could be made faster by using a look up table approach. Using the 
look up table values as initial guesses, the accuracy in the computed values could be 
improved by Newton's method as 




ar," 


2(7;"), 


where. 


1 N N 


dZjT,) 

K V .-=1 i=. 

X {TRj exp(-ri?j. / T, )(1 - exp(-r/?,. / T, ) [S, (1 - exp(-rR ^ / 7; )) - 5/1 - exp(-rR,. / Tj ))] 

+ exp(-ri?^. IT, )(1 - txpi-TR, IT, ))[SjTR, exp(-rR,. IT,)- SiTRj exp(-TR^. IT,)]}. 

With this procedure, the same set of data wjtfi under same conditions, the computation 
of Tl map has taken 2.31 sec only. 




(e) Ti map 


Figure 2.3: First four are spin echo images with TR/TE 
image is the computed Ti map computed. 


= 850, 1100, 1650, 2200/20 ms. The last 







2.2.2. A Least Squares Method for T 2 Calculation 


In this subsection, a non-linear least squares method for computation of T 2 values 
proposed by Xing Li and J.P. Hornak (1994) is described. The apparent T 2 is usually 
measured using the magnetic resonance signal from a spin-echo sequence, 

S{TEi ) = k exp(-r£,. /T^). (2.2.2) 

The proportionality constant k is dependent on such factors as the spin density, amplifier 
gains, spin-lattice relaxation time Ti, and the repetition time of the imaging sequence. 
The apparent T 2 may be calculated from the above equation with an appropriate set of n 
images recorded with different TEj values. The image represents the actual magnetic 
resonance signal. Si, at time TEi , which is composed of a contribution from the spins, 
S(TEi), and random additive noise N, independent of the S(TEi). 

Si =S(TEi) + Ni. (2.2.3) 

Xing Li and J.P. Hornak (1994) studied two techniques for finding the best S(TEi) 
match with the data based on a linear and a rapid nonlinear least-squares solution to 
equation (2.2.2). 

The linear procedure linearizes equation (2.2.2) by taking the logarithms of both 

sides 


ln(Si) = lnik)-TEi/T, 

such that the slope of a line obtained from plotting In(Si) versus TEi is -I/T2. The 
procedure seeks to minimize (t)(T 2 ) by setting 


: 0 and 


30 


30 

— \j aiiu ^ 

dT, 3(ln(S,.)) 


= 0 


for the n data points, where 

1=1 

The slope obtained from a linear least-squares solution is 


?=1 z=l 


i=l 





The difficulty with this method lies in the minimization of <|). The deviation of In(Si) 
from ln[S(TEi)] at each i is treated with equal significance. Taking logarithms of both 
sides of equation (2.2.3) yields 


ln(5,)-ln[5(r£,)] = ln 


1 + 


S{TE,) 


The above equation shows that ln(Si)-ln[5(7’£',)] is dependent on the S/N, which itself 
dependents on TEi. The presence of small signals in the data set used to calculate T 2 
with a linear least squares procedure could severely affect the accuracy of the fitting. 

Alternatively, T 2 could be estimated with a nonlinear regression by minimizing 


HT„k) = ^{S,-kcxp(-TE,/T^)f . 

i=i 


Therefore 

( 50 ( 7 ,, k) = ^5k+^5T2=0. 
^ dk dT^ 


Because neither dk nor 6 T 2 may be zero, then 


^♦=0andi^- = 0. 


dk 


dT, 


These equations give 

2(5, • -koxp(-TE, /T 2 )Qxpi-TE, /T^) = 0, 

1=1 


and exp(-7E,. / ) }7B,. exp(-T£,. / 7^ ) = 0 . 

(=1 

Solving these equations for k and equating both equations yields 

YS, exp(-7£, /T 2 ) J^S^TE, exp(-7E, /T^) 

Jexp(-27E,. /T 2 ) ^TE, exp(-277,. IT^) 

,•=1 i=l 

Now this equation is a function of T 2 only. After reordering the summations one obtains 
a function Z(T 2 ) given as 

exp(-7£, IT 2 )TE^ exp(-27£^. IT^) 

i=l ;=1 

- SjTEj exp(-7E,. /72)exp(-277, IT^)]. 

Xing Li and J.P. Hornak (1994) have solved this equation using a controlled iterative 



calculations is less sensitive to the presence noise in the images in the data set compared 
to linearized least squares method. Since this method requires images acquired at 
different TE's with same TR values, one could use multi echo sequence images without 
any excess scanning time. 

Figure 2.4(a)-(d) shows the multi spin echo images of a volunteer acquired with 
TE/TR = 20, 46, 80, 1 14/2200 ms. Figure 2.4(e) is the T 2 map computed from Figure 
2.4(a)-(d) images using the above method. White dots outside the brain T 2 map are due 
the presence of noise more than the threshold value, here it is taken as 60, in the input 
images. Computing T 2 map with this data using the method described above has taken 
37.31 sec on a Pentium II P.C. with 350 MHz speed and 96 MB RAM. 

This method could be made faster by computing T 2 map using two point method 
with two images from the data set (preferably one image with minimum TE and other 
one with maximum TE to have less errors due to the presence of noise). Now using 
these values as initial guesses, the accuracy in the computed values could be improved 
by Newton’s method as 


rjptl^X 


dT2 


Z{T^), 


where. 


exp(-TE,. /T2)exp(-2TE^. IT^){TE, +2TEj) 

^^2 T 2 i=i j=i 


- Sj expi-TEj /r2)exp(-2TE,. /T 2 )iTEj + ITE ^)} . 

With this procedure, the same set of data under same conditions, the computation of T 2 
map has taken 12.09 sec only. 




(e)T2map 

Figure 2.4: First four are multi spin echo images with TE/TR = 20, 46, 80, 144/2200 ms. The last 
image is the computed Ta map 
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37, 

37, 


37, 

372 

3p 

372 

37J 

372 

37, 

37, 

3p 

373 


37, 

37, 

37, 

3p 


= ^^exp(-(r£i +TE2+TE^)/T2)txp{-(TR^ +TR 2 +TR2)/Ti) 

T'l 7^2 

TR,-(2TR, +l-2/^''2''>) 
xTR2-(2TR2-TE2)e'^^^-'^'^^ TE^ie'^^^''^' +\-2e'^^^'^'^^) +i-2/^^/27-, ^ 

r/?3 -(2r/?3 +1-2^''^^ '''’■) +l-2e^^3/2r,) 

=-^exp(-(r£i +r£2 +r£3)/r2)exp(-(r/?] +r/?2 +7’7?3)/7i) 

T{T{ 

X { [77?, - (27/?, - re, )e^^> WE^ - TE^ )a2a3 
+ [7/?2 -{2TR2-TE2)e’^‘^-'^'^'WE2 -7E,)a,a3 
+[r/?3 -(2r/?3 -7£3)e”'''^^'](r£:, -r72)«ifl2} 
where +1-2/^'^^"'. 

Therefore 

I A 1?^ 0 ^ yt,r/?, (TBj - 7^3 ) + k^TR^ (TE^ - 77, ) + ^37^3 (77, -TE^)^^, 
where it, =[l-(2-77, /77,)/^'^^^‘]a2a3, 
k2 =[l-(2 - 772/772)/^^^^^‘]a,a3, 
it 3 = [1 - (2 - 773 7773 )e^^ 3 / 2 r, 

Therefore, except for some physical situations where the determinant of A vanishes, it is 
possible to compute Ti, T2, and p simultaneously. 


2.2.4. Reduction to a Single Variable Case 
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Most of the methods for tissue parameters T., T, and p estimation in MR imaging either 
require very specific data (e.g., in two point fit (J.P. Jones (1988b) for Ti -calculation, 
where one uses same TE but one of the TR is double of the other and in the iterative 
method (J. Gong and J.P. Hornak (1992)) for Ti -calculation data at different TR's which 
takes more scanning time is required) or are based on solving a possibly over- 
determined system of nonlinear equations in three variables (H. Iwaoka, T. Hirata and 
H. Matsuura (1987)). In the former case it is only approximated signal expressions 
which are used and so some corrective procedures are still to be used. The algorithms in 
the latter category suffer from the well known problems arising in the case of several 
variables, like the close choice for the initial guess and the convergence problems in 
genet dl aiising due to the high degree of freedom available in dimensions larger than 
one. However the situation in the case of a single variable is relatively simple, where the 
initial guesses as well as monitoring the convergence could be affected much more 
conveniently. In this subsection we introduce a method based on separating the 
variables so that finally we end up with an equation of a single variable Ti, whose 
determination subsequently leads to finding the other parameters Ta and p. 

The magnitude signal of spin echo can be expressed as 

Sj = kp(l-2expi~{TRj -TEj /2)/T,) + exp(-rRy ITy))QX^{-TEIT 2 ) 

If we have three independent images, considering the following ratios, 

f = cxpHre, -TE, )IT,)X (l■^^xp(-ra,/^,)-2exp(-(re, -r£,/2)) 

S 2 " (l + exp(-77?2/7’i)-2exp(-(rR2-71E2/2)) 

and 


*^3 

^2 


exp(-(7’E^ 

^ ^ (l + exp(-rR2/7;)-2exp(-(77?2-7^£2/2)) 


we can write 






(.te.-te,) 




S 


2 J 


(1 + exp(-rRi ITi ) - 2exp(-(ri?i - TE^ 

(1 + exp(-77?3 / Ti ) - 2exp(-(rR3 - TE^ / 


X (1 + exp(-rR2 ITi ) - 2exp(-(rR2 - TE^ 


Now, 
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^ |i[(i + exp(-r/?./ri)-2expHr/? -TE^ , 

;=l ;=1 

when 


j = 1, ;; = 2 , 7 ;^ = 3; 

2 “ 2, 7 ^ — 3, 7 ’^^ = 1; 

j ~ j+ ~ I 5 y++ ~ '^• 


Taking logarithms on both sides yields 


F{T,) = Y(TEj,-TEj,,)\n 
;=i 


1 + expi-TRj / r, ) - 2 exp(-(ri?^. -TEj 1 2) I T^) 


Sj 


= 0 , 


and 


J_y iTEj+ -TEj^J 

ar, " 7f 1 + expi-TRj / r, ) - 2 txpi-iTRj - TEj / 2) / T, ) 


xiTRj exp(-TRj /Ti)-2exp{-(TRj -TEj/2)fT^)iTRj -TEjl2)) . 

Starting with the initial approximation 7,°, iteratively Ti values can be computed using 
Newton’s method as: 


^n+l 


= 7 ;" - 


^ dF 


ar," 


F(T"). 


In case of close initial guess problem, a Look-up Table can be used. 

It is observed that at all routinely used TE, TR’s and over the range of clinically 
interested Ti values, the function F(Ti) is, numerically, monotone and guarantees the 



Graph 2.3: Plot of F(Ti) vs Ti values with TR/TE = 2200/20, 1024/14, 2200/80 ms 
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non vanishing derivative and consequently the above iterative method is meaningful. As 
an illustration, Graph 2.3 shows the F(Ti) against the T, values for the spin-echo 
images with TR/TE = 2200/20, 1024/14, 2200/80 ms. 

Now, T 2 and p can be computed as follows: 

Considering two images at a time and substituting the calculated Ti values in the signal 
intensity expression we can have calculated T 2 values for each pixel as 



Note that this expression is different from the conventional two point fit method (J.P. 
Jones (1988b)) which uses the assumption that TR»TE and leads to an approximation 
to actual signal expression. The average of all these values can be considered as T 2 
image and it minimizes the sum of the squares of the errors (SSE). And using these Ti, 
and T 2 values from intensity expression, proton density can be obtained from 

p ^ and their average would give p map minimizing SSE. 

8 jiT.JfjiT,) 

This formulation holds for any pulse sequence whose signal intensity can be 
written as: 

First column of Figure 2.5 shows the Ti, T 2 , and PD weighted images. Ti 
weighted image is obtained from single spin echo sequence with TR/TE = 1014/12 ms 
while PD and T 2 weighted images are the first and second echo images respectively of a 
double spin echo sequence with TR/TE = 2200/20,80 ms. The corresponding Ti, T 2 , and 
PD maps, computed using the above method, are shown in the second column of Figure 
2.5. Figure 2.6 shows the divergent map where the iterative method fails to converge. 
Most of these divergent points are seen in the brain cranial region possibly because of 
mixed signals in these pixels due to the partial volume effect. 



w I 





Ti weighted 


Ti map 


PD weighted PD map 

Figure 2.5: Images in the first column are Tj, T 2 , and PD weighted spin echo images with TR/TE 
= 1012/14, 2200/80, 2200/20 ms of a patient with perilesional gliosis around the cystic cavity. 
Second column images are the Ti, T 2 , and PD map computed using the method described in 


T 2 weighted 


Tainap 








Figure 2.6: Divergent map 


Parameter values computed with reducing to a single variable method and two point 
methods as described in section 2.1 are shown in Table 2.9. For this, spin echo signal 
intensity values are simulated at different Ti, Ta, and PD values with TR/TE = 1100/20, 
2200/20, and 2200/120 ms. Parameter values computed using single variable case are 
close to the actual values. Computations with two point Ti and T 2 estimations have 
shown significant difference in the computed values in comparison with original 
parameter values. These computed Ti and T 2 values are substituted in one of the image 
intensity expressions to obtain PD values. These errors are due to the approximated 
expressions used in the two point methods. 




PD 

Single variable case 

Two point methods 



value 

T, 

T 2 

PD 

Ti 

T 2 

PD 

500 

70 

2400 

499.99 

70.00 

2439.99 

510.35 



750 

80 

2760 

749.99 

79.99 

2760.00 

768.10 




120 

2840 

999.99 

119.99 

2839.99 

1028.01 

141.78 

2789.96 

^21 




159.99 

3119.99 

1554.54 

185.82 

3142.51 




1999.99 





3485.49 

3000 



2999.99 

279.99 

3999.99 



4145.43 


Table 2.9: Parameters computed with single variable and two point methods 


If n(>3) acquired MR images with different pulse sequences or sequence 
parameters are available, then it is possible to make "C 3 data triples from which Ti can 
be computed. By considering suitable consistence data triples, a least squares 
formulation can be made as: 



































E 

^ triples 


;^ft),(ln(/,(r,))-ln(5,)) 


which reduces to solving 


S 

triples 




Z=1 


dT, 


;j 


= 0 , 


2.2.5. Least squares method for Ti, I 2 , and p 


Under the ideal conditions, MRI signal expression can be written as 

S{TE,JR^) = pf^{T,)g,{T^). 

Let Si represent the signal observed at {TEi, TR). 

Consider 


(/»(r„ , p) = X iS{TE ^ , TR , ) - 5,. ]^ 

(=i 

where n(>3) is the number of input images. 

Now the minimization of above equation with respect to Ti, T 2 , and p gives: 


d(p 


0 implies 5^[5(r£, .]//(!, )g,.(r 2 ) = 0 = Fi(ri,T 2 ,p), say. 

i~\ 


d<l> 


= 0 implies J^[S{TE,,TR,)-S,]TEJ^{T,)g,{T^)=0 = F^{T„T^,p), say. 

Z=1 


1^ = 0 implies 2[5(TF, .,r/?,)-5,]/,.(ri)g,.(T2) = 0 = F 3 (ri, 72, p), say. 

i=l 

By Taylor series expansion of Fi(i=l,2,3) around some initial approximation Xo(X = 
Xo+M) and neglecting 2"'* and higher order terms, 

AAX = -F;,^, 



r A 

'^arpxo 


V3p-'Xo 





where, A= 

('■^1 


3p ^Xo 

, AX = 

<1 

, and F = 

F2 



*^372^X0 

^ dp j 


[Ap^ 


1^3 J 


For initial approximations, one could use the method described in subsection 2.2.4 by 
considering any three consistent images. This method is different from the one 
mentioned in the subsection 2.2.3 which linearizes the system of equations first and then 
solves the overdetermined system (in the case of n > 3) by a least squares method. 



2.2.6. Computables of Tissue Parameters and Field 
Susceptibility Maps 


Signal intensity expression for most of the pulse sequences at a pixel can be expressed 
as 

kpf{T,)tx^{-TElT^) = S, 

here we assume that the k is slowly varying and may arise due to gradient and signal 
receiver coil imperfections. 

Taking logarithms on both sides yields 

TE 

ln/fc + lnp + ln/(ri)-— = ln5. (2.2.5) 

^2 


Differentiating in the x direction gives 

In 


TF 

\~ I J rp \ ^ ^ 'V ^ 

-p -j 1 -j r-i2;c ~ x'> 

P " /(T.) Ti 


1 

S 


where Sx is the rate of change in the signal, and Px, Tix, T 2 x are rate of changes in proton 
density, Ti, and T 2 respectively in x-direction and denotes 3/(Ti)/3Ti. 

Sx can be calculated either by central divided differences or by direct Fourier 
differentiation of the image. From the well known result in Fourier analysis we have 

= |exp(-/(^)/'(x)dx = -J (-/i^)exp(-ii^)/(-^)^ = f(0 ■ 

Using this the derivative Sx could be computed. 

In the case of divided differences. 


= 


h 


h being the pixel spacing. The values at half pixel width can be obtained by doubling the 
image through zero padding. 

Suppose we have four images of a cross-section with different TE's and TR s, we can 
write 


1 . .IMA 

p 


TE. 




1. 

Si 


5^, i=l, 2,3,4. 


Subtracting first equation from second and third, we have 


,{TE^-T^r. 

i 1 y. n ^ 2x c 

fiiTi) MTO S^ 


^Ix 

Si 


and 
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' m, ) /,'(r,) 




(re, -re,) 


T2,= 




From the above equation, 


Im.= 

1 



t2 

^2 

{TE^-TE,)' 

1^3 

J 


' m,) f;{To 
,MT,) f,(T,) 



Substituting this in the previous equation gives 


ff/aCT’i) 

1 

<N 

1 

f/A) /,'(r,)'|l 

/'(’'■a 

1 

[a ( 7-,) /,(r,)Jj 


( s 

^2x 


TE 2 - 

-TEy 



. “^2 

J 

TE,- 

-TEy 

1^3 

J 


This can be written as 


( 2 . 2 . 6 ) 


8iTi)Tu=D, 


where 


giro 


f2(Tx) m,) 


TE 2 -TEi 
TE^ -TEy 


\r r' 


//^i) fi'iTQ 

fsiTx) MTr) 


>Tu, 


and 


D = 


So. S 


zx 

' 


TE 2 -TE^ 
TE-i -TE^ 


V 


So. s 


A *^3 


3x __ 

5, 


Similarly with another combination of images we can have 

g(T0T,,=D, 


where 


and 


g(X,) = 


{(fz^Tx) /;(r,)] 


(TE2-TE,'Xf;{Ty) /4'(Ti)| 

Ia^,) hVo] 


[re3-re,(/3(r,) /,(r,)Jj 


( s 

^2x 

c ^ 
_ ^4x 


(TE 2 

-TE,] 

( s 

^3x 

C ^ 
^4x 

, Sz 

■^4 ; 


[te^ 

-TE,'j 

1^3 

^4 J 


From these two equations we can have 


8(Tr) _ D 

m) D 

which can be used for the generation of look-up table for Ti values. Using these values 
as initial approximations, the accuracy in the computed values could be improved 
iteratively by Newton's method. 
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Once Ti is obtained, Tix can be computed from any one of these two equations as 

D 


Tu = 


giro 


Substituting these Ti and Tix values in the equation (2.2.6), T 2 x/T 2 ^ can be computed and 
using these values, the relative proton density variation, px/p, can be obtained. In the 
following system of equations formulation, we show that just increasing the number of 
independent images would not help in computing T 2 , T 2 x, p, Px, and proportionality 
constants. 

Consider three consecutive pixels at i-1, i, i+l in four given independent images. 
Without loss of generality we can take ^i=l and the equation (2.2.5) for four images in 
this case can be written as 

r r.i \ 

Inp'^ 



In k. + In p 

where p^ ,T^ , f ^ {T^), and represent the proton density, T 2 relaxation time, 

function of Ti, and acquired signal at the location j in all four images and j = i-l, i, i+l. 
These equations can be written in the matrix- vector notation as follows: 


0 

0 

0 

1 

-TEi 

0 

0 

0 

0 

1 

0 

0 

1 

-TE^ 

0 

0 

0 

0 

0 

1 

0 

1 

-TE, 

0 

0 

0 

0 

0 

0 

1 

1 

-TE^ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-TEi 

0 

0 

1 

0 

0 

0 

0 

1 

-TE^ 

0 

0 

0 

1 

0 

0 

0 

1 

-TE, 

0 

0 

0 

0 

1 

0 

0 

1 

-TE^ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-TE^ 

1 

0 

0 

0 

0 

0 

0 

1 

1 
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to 

0 

1 

0 

0 

0 

0 

0 

1 

-TE-^ 

0 

0 

1 

0 

0 

0 

0 

1 

-TE. 


ln^2 ^ 

In 

In k^ 
Inp'^ 
l/7-« 


Inp' 

\ItI 



"infcr7/r‘(r,)^ 

insr'//r(j"i 

insrV/r'cT'i 

in(sr//r'(7’.) 

in(s;//i(r,) 

infc^7/“(r,| 

lnSf7/3'«(r,) 

,in(sf7/77r,( 



By the following elementary row and column operations one can see that the column 
rank of the above overdetermined system is less than nine which is the number of 
unknowns to be determined. 


R2-R1 

Re - R5 

Rio - R9 

C5 + TE1C4 

C5 + (TE2-TE0C, 

R3-R, 

Ry - R5 

Rii - R9 

C7 + TEiCe 

C 5 + (TE 3 -TE,)C 2 

R4 - Ri 

00 

1 

Ri 2 - R9 

C9 + TEiCg 

C5 + (TE4 - TEi)C3 


C5 + C7 + C9 = 0. 

Therefore, 

(i) the rank is not nine, 

(ii) if any two of TEi, TE2, TE3, TE4 are distinct, then the rank is eight, and 

(iii) if TE's are all equal, the rank is six. 


Suppose n images are given, then the coefficient matrix size will be 

3n X ((n-l)+6) = 3nx{n + 5). 

Applying the following elementary row and column operations on the quotient matrix, 


to 

1 

Rn+2 - Rn+l 

^2n+2 ~ ^2n+l 

Cn+l + TEiCn 

Cn+l + (TE 2 - TEl)Ci 

R 3 - Ri 

B.n+3 ~ R-n+l 

^2n+3 *” ^2n+l 

Cn+3 + TEiCn+2 

Cn+l + (TE 3 - TEi)C 2 

Rn - Rl 

Cn+1 

Rln ■ Rn+I R3n- " R2n+1 

+ Cn+3 + Cn+5 = 0. 

Cn+5 TEiCn+4 

Cn+l + (TEn - TEi)Cn-l 


So, rank is less than (n+5), which is the number of variables to be determined. Thus, 
just by increasing number of input images and with no conditions on TE's, TR's would 
not help in computing Ti, p, and the factors k 2 ,..., K. 

However, we assume that two of these images have same proportionality 
constants. This assumption is meaningful because, two images from a double echo 
sequence which are obtained under the same experimental conditions except for 
different echo times may have same proportionality constants. By taking ratio of these 
two, T 2 can be computed as 


[StfxiT,)] 

Ti values used here are computed from above formulation which are free from the effect 
of proportionality constants. Substituting Ti and T 2 in the intensity expression in which 
proportionality constant assumed as one, p can be computed. Then all other parameters 



DO 


If whole image is multiplied with a single constant (say, k^) the following least 
squares estimation could be used for evaluating the constant after computing Ti, T2 and 
p as above. 

N 

mmY,[k 2 -S^if =(l> (say), 

^3 i=i 

where 5^3 = /(f 2 ,iT^)exp(~TE 2 /T 2 )p) , and N is the total number of pixels in the 

dd) ^ ~ 1 ^ ~ 

image. — — = 0 implies ^ [^3 - 53,- ] = 0 . That is ^3 = — ^ >?3; , which is the average of 
^^3 ;=i N 

the constants computed at each pixel. 

For simulation experiments, four spin echo images are synthesized with TR/TE 
's as 2200/15, 1550/20, 4650/12, and 2200/120 ms. Two images are taken with same TR 
and different TE's. Table 2.10 shows the % relative errors in the computed Ti, T2, and p 
values due to the approximations when divided differences are used for evaluating 
derivatives. Figure 2.7 shows the maps of Tu, T2x/T2^, Px/p and their histogram 
equalized maps. 



Ti values 

T 2 values 

PD values 

% Li error 

5.586293 

0.091141 

1.895488 

% L 2 error 

16.126315 

5.224453 

9.693125 

% Loo error 

86.1 

50.0 

86.525781 


Table 2.10 : % Errors in computed Ti, T 2 , and PD values when derivatives are evaluated by 
divided differences. 

One of the images is multiplied with a constant 1.1 and parameter values are 
computed using the method discussed in subsection 2.2.4. Figure 2.8 shows the 

parameter maps computed without and with constant multiplication. Errors in the 

computed Ti, T2, and p values are tabulated in Table 2.11. But, computation of 
parameters with the present method did not show any differences in the parameters 

when computed with and without constant multiplication. This is clear from the 

formulation. 


K = l.l 

Ti values 

T 2 values 

PD values 

% Li error 

29.9369 

10.4630 

16.8470 

% L 2 error 

26.6413 

33.4192 

17.0726 

% Loo error 

25.3749 

85.8198 

34.0824 


Table 2. 11: % Errors in computed Ti, T 2 , and PD values with one.of the images multiplied with 




To create small variations in images, each image is multiplied with varying 
factors. These factors are generated from the following computational procedure with 
different k values. 

do i = 0 to 

tempi = A: O' - M2) W 
do 7 = 0 to 

temp2 = k(j-N Ilf/ 
factor = 1 + l/(l+templ+temp 2 ) 
end do 
end do 

For the present study, we have taken A^ = 5 12. For different values of k, one can generate 
different factor maps. 

Figure 2.9 shows the maps of proportionality constants computed as above with 
^ = 1, 2, 3, and their histogram equalizations. As these factors are small, each input 
image is multiplied with (1+factor). Tissue parameters are computed in the presence of 
these factors with ky =^ 14 = 1 , and ^2 = 2 & ^ 3 = 3 . Errors in the computed parameters are 
shown in Table 2.12. 



Ti values 

T 2 values 

PD values 

% Li error 

0.084175 

0.018334 

0.192346 

% L 2 error 

0.515941 

0.050305 

0.223677 

% Loo error 

26.13333 

0.328222 

8.951550 


Table 2.12 : % Errors in computed Ti, T 2 , and PD values using images multiplied with factor 
for different k values 


Table 2.13 shows the errors in the estimation of varying multiplication factor ( for ^2 & 
ks). These high errors due to the approximated derivative evaluation and also, these 
multiplication factors are not locally constant which the present formulation assumes. 



Factor for k =2 

Factor for k = 3 

% Li error 

39.985179 

47.348269 

% L 2 error 

58.868919 

62.736832 

% Loo error 

79.385237 

84.286345 


Table 2.13: : % Errors in estimation of multiplication factors 
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Figure 2.7: Maps of Tjx, T2x/T2^, Px/p and their histogram equalizations 









Ti map without multiplication factor Ti map with multiplication factor 


T 2 map without multiplication factor 


T 2 map with multiplication factor 


PD map without multiplication factor PD map with multiplication factor 

Figure2.8: l®' column: Computed parameters from three weighted images. 2"'’ column: 
Computed parameters from the same images except for one image is multiplied with 1.1 
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Figure 2.9: Different factors and their histogram equalized maps used as simulated varying field maps 





2.3. Error Propagation in Tissue Parameters Estimation 

MRI signal expression can be written as 

The error propagation through these equations into the estimated T\, T 2 , and p can be 
found from the total differential of function ^ of three independent variables Ti, T 2 , and 
p and when applied to three sequences gives 


dp 


»■ -t*-' -t*' 

The errors in 0/ affects the computed Ti, T 2 , and p. These can be computed using 
Cramer's rule as 

(5Ti = — , 5T^= — , 6p= — , 

^ D ^ D D 
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In most of the pulse sequences g{T 2 ) takes the form expC-TE'/Tz). Using this, the 
errors in estimated Ti, T 2 , and p values can be expressed as 
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2.3.1. Errors due to Bit Reduction 

The sensitivity of the tissue parameters computations are checked numerically by 
calculating tissue parameters from images with different depths. The depth (number of 
bits used to store image information) is reduced at each level by one bit. Images in the 
Figures 2.10-2.12 are obtained from Ti, Ta, and PD weighted spin echo images of a 
normal sample by reducing bits successively from them. Visually these images do not 
show any significant changes up to five bits reduction. In the seven bit reduction 
deterioration is clear. Figures 2.13-2.15 are the corresponding tissue parameter maps 
computed with the method described in the subsection 2.2.4 using respective images 
from Figures 2.10-2.12. Tables 2.14-2.16 show the errors in weighted images due to 
reduction of bits. The sensitivity of tissue parameters calculation on the depth of images 
used is shown in Tables 2.17-2.19. The computed errors are; 
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The computed average Li, L 2 , and L=o errors in Ti, T 2 , and PD weighted images 
due to reduction of bits is shown in Graphs 2.4(a)-2.4(c) respectively. In all the three 
types of errors, these graphs reveal small amount of errors up to four bits reduction. In 
all three weighted images error magnitudes are almost same up to five bits reduction. 
Five bits onwards PD weighted image has shown more error compared to Ti weighted 
image. T 2 weighted image has slightly less errors compared to PD weighted image but 
has more errors with respect to Ti weighted image. 

Graphs 2.5(a)-2.5(c) show the average Li, L 2 , and L«, errors respectively in Ti, 
T 2 , and PD maps computed from the successively bits reduced weighted images. In all 
the three types of errors T 2 map has shown relatively low errors compared to errors in T 1 
and PD maps. 

Following R.E. Hendrick (1999), the standard average Ti values of 
methomoglobin, white, lower cystic fluid, gray and CSF are 460, 510, 720, 760, 2650 
ms. The corresponding T 2 values are 106, 67, 43, 77, 280. From Tables 2.17 and 2.18, 
the removal of five bits onwards could induce average errors of magnitude around 94 in 
Ti and around 3 in T 2 , which mean that there would be problems in differentiating 
methomoglobin and white matter as well as lower cystic fluid and gray matter. Similarly 
on the T 2 map the white and the gray matter regions could get mixed up. The 
corresponding average L 2 errors were as high as 8.5 and 269.7 for T 2 and Ti 
respectively. Thus a depth of the images below 12 bit which is the current DICOM 
standard would severely affect the use of the parameters in relaxometry studies. Our 
recommendation would be to compute the gray values directly from the machine data 
and not to rely on the machine stored images, if the relaxometry study is to be 
conducted. 
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Figure 2.10: Tj weighted image with different number of bits removed 
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Figure 2.11: T 2 weighted image with different number of bits removed 
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Figure 2.12: PD weighted image with different number of bits removed 
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Figure 2.13: Ti map computed from different number of bits removed weighted images 
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Figure 2.14: Ta map computed from different number of bits removed weighted images 
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Figure 2.15: PD map computed from different number of bits removed weighted images 








No. of bits 
removed 

Average 
Li Error 

Average 
L 2 Error 

Loo Error 


Average 
Relative Li 

Average 
Relative L 2 



0.702713 


■ 

0.018769 

0.045069 

Two bits 

1.499054 

1.867604 

3.000000 


0.062837 

0.147732 



4.177157 

7.000000 


0.139810 

0.271089 

Four bits 

7.488068 

8.768975 

15.000000 


0.275875 

0.452932 

Five bits 

15.498810 

17.843015 

31.000000 


0.484417 

0.667783 

Six bits 

26.406036 

30.854564 

63.000000 


0.610945 

0.757743 

Seven bits 

40.733185 

51.932167 

127.00000 


0.667506 

0.783768 


Table 2.14: Average & average relative errors due to bit reduction in Ti weighted image 


No. of bits 
removed 

Average 

Li 

Average 

L 2 


1 

Average 
Relative Li 

Average 
Relative L 2 

One bit 

0.499908 

0.707042 

1.000000 


0.007939 

0.015571 

Two bits 

1.493866 

1.866149 

3.000000 


0.023826 

0.041362 

Three bits 

3.482758 

4.166777 

7.000000 


0.055594 

0.092985 

Four bits 

7.463226 



■ 

0.119498 

0.195517 

Five bits 

15.522552 



■ 

0.250741 
1 

0.394480 

Six bits 

33.835052 

38.111704 

63.000000 

■ 



Seven bits 

54.470795 

61.542353 

127.00000 


0.665562 

0.786245 


Table 2.15: Average & average relative errors due to bit reduction in T2 weighted image 



BSSIBBESBi 

0.706686 

1.000000 

0.006696 

0.013916 


1.876468 


0.020167 

0.036777 

Three bits 

3.497849 

4.179638 

7.000000 


0.082950 

Four bits 

7.494675 


0.100701 

0.174544 

Five bits 

15.468552 


0.210812 

0.355055 

Six bits 

33.158981 

37.763136 

63.000000 


0.443592 

0.635873 

Seven bits 


0.591845 

0.748702 


Table 2.16: Average & average relative errors due to bit reduction in PD weighted image 






























































































Table 2.17: Average & average relative errors due to bit reduction in Ti map 



No. of bits 

Average 

Average 

IIBSHI 

Average 

Average 

removed 

Li 

U 

innnni 

Relative Li 

Relative L 2 


0.001856 


0.003125 


Two bits 

0.398558 

1.433227 

90.617401 


0.004496 

0.007441 

Three bits 

0.812391 

2.629004 

200.93664 


0.009430 

0.014707 


4.958928 

231.10553 


0.019176 

0.029601 

Five bits 

3.227306 

8.553119 

352.02179 

■—a 

0.055804 

Six bits 

5.981129 

11.752798 

349.32559 


0.073766 

0.099153 

Seven bits 

10.592555 

19.709234 

399.02548 


0.123114 

0.158454 



Table 2.18: Average & average relative errors due to bit reduction in T 2 map 



No. of bits 
removed 

Average Li 

One bit 

3.238593 

Two bits 

8.008484 

Three bits 

17.757535 

Four bits 


Five bits 


Six bits 


Seven bits 

231.608932 



6612.5727 


Table 2.19: Average & average relative errors due to bit reduction in PD map 


0.001781 

0.006280 

0.004402 

0.017449 

0.009669 

0.040801 


0.092453 


0.091865 

0.283429 

0.161765 

0.366470 
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Graph 2.4(a): Average Li errors in Ti, T 2 , and PD weighted images due to bit reduction 



Graph 2.4(h): Average L2 errors in Ti, T2, and PD weighted images due to bit reduction 



Graph 2.4(c): Average errors in Ti, T 2 , and PD weighted images due to bit reduction 
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'aph 2.5(a): Average Li errors in Ti, T 2 , and PD maps computed from bit reduced images 



raph 2.5(b): Average L 2 errors in Ti, T 2 , and PD maps computed from bit reduced images 



Graph 2.5(c): Average L„ errors n Ti, T 2 , and PD maps computed from bit reduced images 


2.3.2. Effect of Filtering 


All MR images slightly, sometimes severely, suffer with filtering (blurring) effects due 
to the image acquisition techniques and reconstruction methods. It is also a common 
practice in routine clinical diagnosis to apply some filters to minimize the noise in MR 
images. Due to the application of filters, signal also gets affected and if these filtered 
images are used for tissue parameters calculation, the parameter maps would get 
affected. To demonstrate this numerically, we study the effect of Poisson filter in the 
calculation of Ti, T 2 and PD maps. 

Poisson sums are obtained by convolving / with the Poisson kernel. For 0 < r < 1, the 
Poisson kernel is given by 

“ . 1 — 

P{r, r) = 1 + 2 V r-^ cos jt = 

j=\ l-2rcosr + r^ 

(/) = — + ^ r* (a^, cos kx-vbp. sin kx) . 

k=i 

Here filter coefficients are generated by using (0 < r <1). In the present numerical 
studies we have taken r = 0.99. 

Recently, R.K.S. Rathore, R.K. Gupta, V.S.N. Kaliprasad, S.B. Rao, and 
Sushmita Datta (2000) have proposed a de-blurring technique to improve image 
resolution by removing or minimizing the effect of filters. Images used after resolution 
enhancement for parameter maps calculation gave better contrast maps. For study of 
errors, first we have computed Ti, T 2 , p values using three weighted images of a 
volunteer. Using these values three spin echo weighted images with TR/TE = 1012/14, 
2200/20,80 ms are synthesized. These images are considered as original images and 
associated tissue parameter values as original Ti, T 2 and PD values. Synthesized 
weighted images are blurred and iteratively de-blurred by Poisson sums. Figures 2.16- 
2.18 show the blurring effect and de-blurring of the three weighted images by Poisson 
filter. In these three Figures first image is the actual synthesized weighted image and 
second one is obtained by filtering the synthesized image by Poisson filter. Third, 
fourth, fifth, and sixth images are the de-blurred images corresponding to first, fourth, 
seventh, and tenth iterations of the de-blurring scheme. Visually one could observe that 
these images are gradually approaching to the actual synthesized image as the number of 
iterations increase. Figures 2.19-2.21 are the tissue parameter maps computed at 
different iterations of the de-blurring of weighted images. In all these Figures first image 



represent the tissue parameter map computed from the actual synthesized image and the 
second one from the filtered images. Third, fourth, fifth and sixth images are the 
parameter maps computed from first, fourth, seventh and tenth iterations of the de- 
blurred weighted images. Even in the case of tissue parameters maps, visually, they 
appear to converge to the maps corresponding to the original synthesized images. 

At each iteration Ti, T 2 and PD maps are calculated and errors are computed. 
Errors are computed according to: %error = 100 x || original-computed || / 1| original ||, 
the norm being the Li, L 2 , and Lc« norm. Tables 2.20-2.22 show the errors in weighted 
images due to filtering and iteratively removing its effect. In all three images, in the first 
iteration itself the error is significantly reduced and in the later iterations it is gradually 
reduced. Tables 2.23-2.25 show the errors in tissue parameters computed at each 
iteration of the de-blurring of the weighted images. From these Tables it is clear that the 
errors in computed tissue parameters are less with less filtering effect. So, it is strongly 
recommended not to use any kind of filters while reconstructing the weighted images 
from acquired raw data from MRI scanners if one wishes to compute tissue parameters 
from them. 

Graphs 2.6(a)-2.6(c) show the average Li, L 2 and Lc errors respectively in 
weighted images by blurring and de-blurring by Poisson filter with r = 0.99. From these 
graphs it is clear that as number of iterations increased the errors gradually decreased. 
Graphs 2.7(a)-2.7(c) show the average Li, L 2 and Loc errors respectively in Ti, T 2 and 
PD maps computed from blurred and iteratively de-blurred Ti, T 2 and PD weighted 
images by Poisson filter. Graph 2.7(a) shows decrease in Li errors in parameter maps 
computed from higher iterations of de-blurred weighted images. At all iterations, Ti 
map has shown more Li errors compared to T 2 and PD maps. PD map has relatively low 
Li errors compared to Ti and T 2 maps. Graph 2.7(b) shows the decrease in the L 2 errors. 
T 2 map has shown a significant decrease in error from iteration 5 to iteration 6. In L 2 
errors T 2 map has more errors compared to Ti and PD maps. Graph 2.7(c) shows the L«, 
errors in the computed maps. 
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Figure 2.16: Blurred and iteratively de-bluriing by Poisson filter (Ti weighted image) 
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Figure 2.17: Blurred and iteratively de-blurred by Poisson filer (T 2 weighted image) 











Original 


Blurred 



Iteration 1 


Iteration 4 


Iteration 7 


Iteration 10 


Figure 2.18: Blurring and iteratively de-blurring by Poisson filter (PD weighted image) 
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Figure 2.19: Ti map calculated from weighted images blurred and de-blurred by Poisson filter 
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Figure 2.20: map calculated from weighted images blurred and de-blurred by Poisson filer 
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Figure 2.18: PD maps computed from weighted images blurred and de-blurred by Poisson filter 








Li Error 


L 2 Error 


L„ Error 


Blurred 


15.920006 


47.029857 


95.099185 


Iteration 1 


Iteration 2 


9.584898 


38.124664 


6.788624 


31.447658 


87.723128 


84.033179 


Iteration 3 


5.082074 


26.142247 


80.211153 


Iteration 4 


3.941934 


22.289792 


75.993708 


Iteration 5 


3.155377 


19.297435 


71.236598 


Iteration 6 


2.566577 


16.853122 


65.905023 


Iteration 7 


2.108289 


14.709351 


60.056701 


Iteration 8 


1.765043 


12.837680 


53.827110 


Iteration 9 


1.486457 


11.109766 


47.410174 


Iteration 10 


1.258534 


9.594922 


41.030135 


Table 2.23: % Error in deblurring by Poisson filter (r = 0.99) (Ti map) 


Li Error Lz Error Error 

Blurred 12.611589 52.643458 95.974726 

Iteration 1 7.054603 45.525262 91.927930 

Iteration 2 4.965641 41.441344 86.610244 

Iteration 3 3.841923 37.894183 80.440234 

Iteration 4 3.025411 34.243302 73.746165 

Iterations 2.463599 31.035154 66.638912 

Iteration 6 1.882758 21.886574 64.735094 

Iteration 7 1.507131 18.421377 58.996586 

Iterations 1.288753 17.660786 50.573965 

Iteration 9 1.112232 17.002172 45.170086 

Iteration 10 0.969051 16.410773 44.383822 

Table 2.24: % -Error in deblurring by Poisson filter (r = 0.99) (T 2 map) 


Li Error 


L 2 Error 


Loo Error 


Blurred 


8.434432 


21.572942 


84.377495 


Iteration 1 


Iteration 2 


4.715807 


17.177758 


3.367360 


14.361484 


76.940832 


70.786209 


Iteration 3 


2.562864 


12.052858 


64.917243 


Iteration 4 


Iteration 5 


2.032226 


10.383647 


1.641415 


8.996978 


59.165056 


53.534160 


Iteration 6 


1.345781 


7.831474 


48.075080 


Iteration 7 


1.115034 


6.820065 


42.849856 


Iteration 8 


0.935416 


5.932835 


37.917754 


Iteration 9 


0.792656 


5.146250 


33.327236 


Iteration 10 


0.674427 


4.452812 


29.112363 


Table 2.25: % -Error in deblurring by Poisson filter (r = 0.99) (PD map) 
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Graph 2.6(a): Average Li errors in Ti, T 2 , and PD weighted images blurred and 
iteratively deblurred by Poisson filter (r = 0.99) 



Graph 2.6(b): Average L 2 errors in Ti, T 2 , and PD weighted images blurred and 
iteratively deblurred by Poisson filter (r = 0.99) 



Graph 2.6(c): Average errors in Ti, T 2 , and PD weighted images blurred and 
iteratively deblurred by Poisson filter (r = 0.99) 
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Graph 2.7(a): Average Li errors in Ti, Ta, PD maps computed from blurred and 
iteratively deblurred images 



Graph 2.7(b): Average La errors in Ti, Ta, PD maps computed from blurred and 
iteratively deblurred images 



Graph 2.7(c): Average L„ errors in Ti, Ta, PD maps computed from blurred and 
iteratively deblurred images 


Chapter 3 

Generation of High Contrast Images for 
Diagnosis 


For medical diagnosis adequate contrast among normal tissues is necessary for good 
anatomical definition, and good contrast between normal and diseased tissues is 
essential for sensitivity to disease. The primary sources of inherent tissue contrast in 
MRI are threefold - hydrogen spin density, Ti, and T 2 relaxation times. Although 
hydrogen densities within soft tissues typically vary by a few percent, the proton density 
contribution to the measured MR signal tends to vary by a greater amount, up to 30% 
among soft tissues. T 1 and T 2 relaxation times often vary even more widely, sometimes 
varying among soft tissues by more than 100%. For example, the Ti and T 2 values of 
CSF are several times the Ti and T 2 values, respectively, of white matter and gray 
matter (R.E. Hendrick (1999)). 

Usually, Ti weighted images show good anatomical details while T 2 weighted 
images show a good contrast between white matter, gray matter, and CSF. T 2 weighted 
images are generally used for pathological observations. Even though PD weighted 
images are relatively flat, some pathologies like, perilesional gliosis around cystic cavity 
are better visible on PD weighted images than on Ti and T 2 weighted images. 

This chapter mainly deals with the generation of high contrast images. In section 
3.1, an entropy-like maximization procedure which computes the weights is presented. 
Using these weights, a weighted linear combination of the acquired images is produced. 
Given a set of linearly independent images, Gram-Schmidt orthogonalization process is 
discussed in section 3.2 to generate orthonormal images. In section 3.3, generation of 
contrast images through functional maps is presented. Section 3.4 deals with a least 
squares method to produce tissue parameter maps contrast through a weighted 
combination of the images. In section 3.5, generation of good contrast in perilesional 
gliosis region on synthesized IR images is demonstrated. Synthesis of image values at 
divergent points, where tissue parameter values could not be computed due to partial 
volume effect, is discussed in section 3.6. The utility of tissue parameter maps for 
segmentation is shown in section 3.7. 



3.1. Entropy-Like Maximization 


Given p images corresponding to a cross-section, their weighted linear combination can 
be written as 




k=\ 


where Ck's are the weights and S* is the signal intensity in the kth image. These weights 
could be obtained by maximizing the function (j) with respect to Ck's and X defined as 
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where A is a Lagrange multiplier and N is the number of pixels in one image. 
This maximization process gives a system of equations given as 
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The above nonlinear system of equations can be solved iteratively by Newton's method 
as 

C" = C' - A-'F, (3.1.1) 

where 
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Figure 3.1 (a)-(c) show the Ti, Ti, and p weighted images. Figure 3.1 (d) is the weighted 
linear combination of the images (a)-(c). These weights are obtained from equation 
(3.1.1). This image appears similar to PD weighted image, but the better appearance of 
perilesional gliosis (black arrow) than in any acquired weighted spin echo images could 
be observed. 

However, we compute these weights by maximizing the following quantity over 
Cj's as follows: 

N k 

max25;c^.5y,.ln(c^.5'^.i), 
i=l j=l 


subjected to =1. Here, N is the total number of pixels and Sji is the ith pixel 

k=i 

intensity of the y'th image. 

Therefore = 0 , y= 1 ,2, . . . ,^ implies 

dCj 

Cj =txp(-il + (^Sjiln(Sji)/^Sji))), j=l,2,3,...X (3.1.2) 

i=i 1=1 

Figure 3.1 (e) is a weighted linear combination of (a)-(c), with weights being computed 
from equation (3.1.2). A good contrast between gliosis region and surrounding brain 


tissue could be seen. 




PD weighted With weights from equation (3.1.1) 



With weights from equation (3.1.2) 


Figure 3.1: Generation of entropy like contrast images using weighted linear combination of 
acquired images 





3.2. Orthonormal Images 


Given a set of N linearly independent images {Sq,S^,...,S , we can construct an 

orthonormal set { Sg , 5, , . . . , } of images which are linear combination of the original 

set using Gram-Schmidt orthogonalization procedure. This sequential procedure is 
described as follows: 


Set So = ' 


^0 


ii. Define y, as the Si minus the projection of Sj onto Sq : 
yj = S, — < Sj , Sq > Sq 

The image y, is orthogonal to Sq by construction 


iii. Set S, = ■ (i.e. normalizing the result of the preceding step) 

II )'i II 

iv. Define yj as the Si minus the projection of Si onto Sq and Sj : 

-Vi = ‘S'i~ < S2,So > Sq- < S2,Si > Si 

V. Normalize: S, = ^ 

ll)’2ll 

vi. Continue this process until S ^_i has been defined 

Computed results are shown in Figure 3.2. Images in Figure 3.2 (a) and Figure 3.2 (b) 
are the Ti weighted and Ti weighted images. Figure 3.2 (c) is obtained by 
orthonormalization of Ti weighted image with respect to Ti weighted image . This 
shows the Ti weighted image without gray matter and CSF. A clear white matter is 
conspicuous on this image. Figure 3.2 (d) is obtained by orthonormalizing Ti weighted 
image with respect to Ti weighted image. This image shows features complement to 
Figure 3.2 (c). 




(c) Ti weighted with respect to T2 weighted 


Figure 3.2: Generation of orthonormal images 



(b) T2 weighted 



(d) T2 weighted with respect to Ti weighted 


Gram-Schmidt orthogonalization process 







3.3. Functional Maps 


Usually MR weighted images are acquired by selecting specific pulse sequences and 
sequence parameters. For example, spin echo images with short TR and short TE give 
Ti weighted images whereas long TR. short TE and long TR, long TE produce PD 
weighted and T2 weighted images. In these weighted images to some extent the 
contribution from other tissue parameters is also there. 

Once Ti, T2, and proton density maps of a cross-section are computed, then it is 
possible to generate the functional Ti map and functional T2 map by computing the 
values of the functions of Ti and T2 respectively. In the case of spin echo sequence these 
functions will be 

/(Ti ) = 1 - 2 exp(-(rR -TEI2)IT^) + exp(-rR IT ^ ) , 

g(r2) = exp(-rE/r2). 

Figure 3.2 shows Ti, T2 weighted images, Ti, T2 computed maps and the corresponding 
Ti, T2 functional maps. The functional Ti map shows excellent contrast between gray 
and white matter with clear anatomical details and intensities are proportional to that of 
Ti weighted image. Even though T2 functional map shows better contrast compared to 
T2 map, T2 weighted image has good contrast compared to both T2 map and T2 
functional map. 






Ti functional map Tj functional map 

Figure 3.3: Generation of Tj and T2 functional maps 


Ti weighted 


T2 weighted 


Ti map 


T2 map 








3.4. Generation of Parameter Map Like Contrast Images 


Using a weighted combination of weighted acquired images, it is possible to obtain 
images with tissue parameter maps like contrast. Once tissue parameter maps are 
computed, these weights can be obtained by minimizing the sum of the squares of the 
differences between the parameter maps and weighted combination of the weighted 
images. The weights, cj's can be computed as follows. 

-|2 

^jSji =^(say), 

where P, Sj , N, k, and cj are the parameter map, j* acquired image, number of pixels in 
the image, and weight associated with j* image respectively. 

3F n [ k ~ 

Now,- — = 0, p = l,2,-.-,k, gives Pi -^CjS ji S =0, p = l,2,...,k, 

i=l L 2=1 J 

k f N N 

which can be written as: ^ ^ SjjSpi Cj = ^ PiSp^ , p = 1,2 , ...,k . 

j=i ,'=1 J ,-=i 

The weights cj 's can be computed by solving the above system of equations. 

Results are shown in Figures 3.4(a) and 3.4(b) The first four images of Figure 
3.4(a) show Ti, T 2 , p, and MT SE Ti weighted spin echo images respectively and the 
last two images are MTR map and the weighted combination with respect to MTR map. 
In Figure 3.4(b) first column shows the computed Ti, T 2 , and PD maps while the second 
column represents corresponding weighted combinations to generate respective 
parameter contrast images. In all these combinations an excellent contrast similar to 
parameters maps is visible. These weights are computed at five consecutive slices each 
with slice thickness 5mm and inter-slice gap 0.1 mm in a single patient and results are 
tabulated in Tables 3. 1-3.4. These results show that the orders of magnitudes of these 
weights over five slice are comparable. That means the average values of these weights 
could be used for all these slices without computing tissue parameters and weights at 
each slice level. This observation could be extended to computing the weights over the 
slices possessing same anatomical structures and making a data base of them which 
could be used to obtain tissue parameters like contrast without actually computing the 
tissue parameters. Generation of these images would be very fast. For a 256x256 image 
it takes 196608 multiplications and 131072 additions. This data base could be made for 
different pathological cases also. 
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Ti Weighted 



PD Weighted 



MTR Map 



T2 Weighted 




Weighted combination w.r.t MTR map 


Figure 3.4(a): Generation of tissue parameter contrast Images by a linear weighted combination 
of images 









T2 Map 


Weighted combination w.r.t. Timap 



Figure 3.4.(b): Generation of tissue parameter contrast images by a linear weighted combination 
of images 
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Ti Contrast 

La Ave. % Error 

Cl (PD wtd) 

Ca (Ti wtd) 

CaCTawtd) 

Slice 1 

24.556313 

669.596252 

-681.283630 

39.475903 

Slice 2 

25.589626 

611.305359 

-693.740967 

60.294865 

Slice 3 

25.249855 

656.256958 

-758.47882 

67.323715 

Slice 4 

26.137049 

622.984802 

-759.524658 

69.148926 

Slice 5 

25.440390 

653.556030 

-748.149902 

44.035568 


Table 3.1: La average percentage errors in weighted combinations with respect to Ti map and 
corresponding weights 


Ta Contrast 

La Ave. % Error 

Cl (PD wtd) 

Ca (Ti wtd) 

CslTaWtd) 

Slice 1 

21.594936 




Slice 2 

25.580469 

21.918722 

-70.840324 

155.125931 

Slice 3 

28.500360 

17.900223 

-78.899765 

196.744797 

Slice 4 

24.503012 

16.844803 

-62.885311 

169.665085 

Slice 5 

27.511993 

17.655508 

-63.075607 

161.335281 


Table 3.2: La average percentage errors in weighted combinations with respect to Ta map and 
corresponding weights 


PD Contrast 

La Ave. % Error 

Cl (PD wtd) 

Ca (Ti wtd) 

C3(TaWtd) 

Slice 1 

9.513255 

411.917389 

-228.914139 

-56.176613 

Slice 2 

10.261456 

374.833130 

-227.364639 

-51.715515 

Slice 3 

10.445120 

393.376923 

-249.356979 

-55.402931 

Slice 4 

9.478690 

370.599182 

-230.857346 

-57.778591 

Slice 5 

9.533550 

383.339325 

-243.544052 

-60.875698 


Table 3.3: La average percentage errors in weighted combinations with respect to PD map and 


corresponding weights 


MTR Contrast 

La Ave. % Error 

Cl (MTSETiWtd) 



Slice 1 

16.440626 

-199.524506 

209.477615 

6.257372 

Slice 2 

17.437050 

-197.515884 

195.639603 

6.594873 

Slice 3 

18.460693 

-213.683929 

217.655014 

7.248658 

Slice 4 

17.822056 

-203.493179 

211.895462 

7.380937 

Slice 5 

18.662191 

-213.680634 

208.736847 

6.506155 


Table 3.4: La average percentage errors in weighted combinations with respect to MTR map 


and corresponding weights 
























































3.5. Synthetic MRI 


One application of the calculation of standard MRI intrinsic parameters is the synthesis 
ot MR images at desired pulse sequences and sequence parameters to obtain optimum 
contrast. MR image synthesis basically consists of three steps: 

(a) acquisition of standard MR images at different TR and TE values, 

(b) computation of Ti, Ta, and p maps, and 

(c) calculation of a synthetic MR image by substituting the computed Ti, Ta, and p 
values into the equation describing the pulse sequence of interest with desired sequence 
timing parameters. 

S.A. Bobman, S.J. Riederer, J.N. Lee, S.A. Suddarth, H.Z. Wang, and J.R. 
MacFall (1985) have shown that synthetic results can resemble actual spin-echo images 
within about 5% of the measured signal level for imaging of the brain. MR image 
synthesis can serve two main purposes: it can be used as an educational tool and it can 
be used for retrospective contrast optimization in clinical diagnosis work. In the 
following, we demonstrate the usefulness of MR image synthesis in obtaining good 
contrast in ROI compared with neighboring tissues. 

Figure 3.5 first row shows Ti, T 2 , and PD weighted images of patient with 
gliosis around cystic cavity. On Ti weighted image the cystic cavity appears as hypo 
intensity (white arrow), but the gliosis is not visible clearly (black arrow) whereas on T 2 
weighted image it merges with the cystic cavity (black arrow). PD weighted image 
shows a better separation of gliosis from cystic cavity and surrounding brain tissue 
compared to Ti and T 2 weighted images. The second row shows the computed Ti, T 2 
and PD maps respectively. A good contrast in these images can observed. The first 
image in the third row, the first magnetization transfer spin echo Ti weighted (MT SE 
TO image, shows a clearer and bigger gliosis (black arrow) compared to PD weighted 
image. The next two images in the same row are two inversion recovery synthesized 
images with sequence parameters TR/TE/TI = 2200/80/1500, 6000/40/2500 ms 
respectively. The pulse sequence parameters were chosen so that these images show an 
excellent contrast between gliosis and cystic cavity (black arrows). Note that the size of 
the gliosis on these synthesized images matches with MT SE Ti image. Since for a new 
study it may not be possible before hand to figure out the pulse sequences as well as 
their parameters which bring out the desired feature of interest, so long as the pathology 
depepi^s only on the tissue parameters, the best strategy seems to be to determine these 


I I I 

tissue parameters first and then pass to different syntheses. The second image in the 
third row is a clear demonstration of the philosophy expounded in the previous sentence. 



Figure 3.5: row: Ti, T 2 , PD weighted images 

2"“* row: Ti, T 2 , PD maps computed using weighted images 

3'“ row: 1 image: MTSE Ti weighted image 

2"'^ and 3^“ images: Synthesized IRSE images at TR/TE/TI = 2200/80/1500, 
6000/40/2500 ms 



3.6. Partial Volume Correction for Synthesis 

In a large, homogeneous test object, regions for Ti and T 2 estimation are available in 
which the entire voxel is filled with the same material. This is almost never the case in 
the living body. Obvious examples of this can be seen in the brain, where the gray brain 
tissue is a thin layer over the white matter, in intimate contact with it (and also with 
CSF). In such a case, a high proportion of voxels in the general gray matter regions will 
also contain, in variable amounts, white matter and/or CSF (M.A. Foster and A. Haase 
(1996)). Also, since the image reconstruction grid is small comparing with the object 
imaged, there would be more than one type of tissue near the boundaries that have 
contributed to signal in one pixel. It is observed that in these mixture signal cases, the 
iterative method presented in subsection 2.2.4 for tissue parameter estimation leads to 
divergent points (see Figure 2.6). 

Suppose, in an MR image, m different neighbouring tissues are contributing to 
the signal of a pixel with different weights (say, (Oj ; j = 1 ,..., m) and we have m 
independent images corresponding to the same cross-section. It is assumed that 
sufficient independent data is available so that the tissue parameters associated with 
these m different tissues could be computed from it. Now the weights (o/s can be 
computed from the following system of equations 

m 

= S,. , j=l,2,...,m, 

M 

m 

subjected to = 1 • 

;=i 

In the case of spin echo sequence, 

a,j = p^.(l-2exp(-(rR,. -TE, I2)/T,j +exp(-rR,. /r,y))exp(-re, /T^j), 

where Ty , T 2 y, and py are the Ti, T 2 , and p mean values of the jth tissue and TE,-, TR; are 
the echo and repetition times of the ith image. Once these weights are computed, it is 
possible to synthesize the pixels values corresponding to divergent points at different 
pulse sequences or sequence parameters. 

Figure 3.6(a) shows the Ti weighted (TR/TE = 1012/14 ms) image which is one 
of the three weighted images used for tissue parameter estimation. Figure 3.6(b) is the 
divergent map where the parameter values could not be computed. Figure 3.6(c) is the 
synthesized image with sequence parameters same as for Ti weighted image. Image 
values are not computed at divergent points. Figure 3.6(d) is same as (c) except for 
correction at divergent points. Assuming three different types of tissue have contributed 



to cause partial volume at each divergent pixel, the weights a^'s have been computed 
and used them in the synthesis. Resemblance of the synthesized pixels at divergent 
points in the cranial part (white arrows) could be observed. By using more images the 
synthesis could be improved. 



(c) Synthesized without partial volume (d) Synthesized with partial volume correction 

correction 


Figure 3.6 : Demonstration of partial volume correction for image synthesis 







3.7. Tissue Parameter Maps for Segmentation 


Segmentation refers to a subclass of enhancement methods by which a particular 
object, organ, or image characteristic is extracted from the image data for purposes of 
visualization and measurement. Segmentation involves associating a pixel with a 
particular object class based on the local intensity, spatial position, neighboring 
pixels, or prior information about the shape characteristics of the object class. 
Segmentation is a central problem of image analysis because it is a prerequisite for the 
majority of analysis methods including image restoration, shape analysis, motion 
detection, quantitative analysis of pathologies, and volume and area estimation. Even 
though MRI images have excellent tissue contrast compared with images produced by 
other medical imaging modalities, there are no clear boundaries between different 
organs. Therefore, the segmentation methods based on thresholding and boundary 
detection would not give satisfactory results. Some statistical methods use pixel 
intensities and its neighborhood information for segmentation. On the basis of normal 
distribution of pixel intensities in different regions (segments), segmentation can be 
done by maximizing the likelihood estimation function (R.K.S. Rathore, S. Datta, 
R.K. Gupta, S.B. Rao, R. Verma (2001)). All these methods depend upon pixel 
intensity and its neighborhood. It is observed that in different segments tissue 
parameters are normally distributed. In Figure 3.7, the histograms of pixel intensities 
of white matter in Ti, T 2 , and PD weighted images and histograms of corresponding 
tissue parameter maps is shown. The normal distribution of tissue parameters within 
the ROI may be observed. This motivates to use tissue parameter maps for 
segmentation methods like MLE. Since, in most of the diseases, the parameter values 
are significantly different from the normal tissues, these maps may give satisfactory 
segmentations. However, there may be some pulse sequences which could produce 
good contrast in the images than in the parameter maps and consequently better 
segmentation. Table 3.5 shows the contrast between some brain tissues in the case of 
parameters as well as in weighted spin echo images. The weighted images are 
considered with TRITE = 500/15, 2200/80, 20 ms. 

Inherent Tissue Contrast 

Tissue contrast derives primarily from inherent p, Ti, and T 2 differences. Other 
sources of tissue contrast are flow, magnetic susceptibility inhomogeneities, and 



chemical shift. Inherent p, Ti, and T 2 contrast between two tissues (A and B) can be 
defined in terms of the characteristics of the two tissues as follows (R.E. Hendrick 
(1999)): 
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_ Pb Pa 
Pb+Pa ’ 
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C. = 


T2 -Ti 


Inherent tissue properties remain essentially unchanged regardless of the imaging 
method. Although defined only for MRI, Ti, T 2 , and p will not change simply because 
a different pulse sequence or a different set of timing parameters is chosen. 



Parameter Contrast 

Signal Contrast (SE) 

Ti 

T 2 

PD 

Ti wtd. 
(500/15) 

T 2 wtd. 
(2200/80) 

PD wtd. 
(2200/20) 

Gray-White 

0.196850 

0.069444 

0.061538 



0.058711 

-CSF 

0.554252 

0.568627 

0.183432 

0.256499 

0.287247 

0.019660 

-Edema 

0.084337 

0.241379 

0.109677 

0.085642 

0.285852 


White-CSF 

0.677215 

0.613833 



0.389265 

0.078281 

-Edema 

0.276596 

0.305699 

0.170068 

0.030490 

0.387974 

0.199471 

CSF-Edema 

0.492958 

0.379310 

0.075269 

0.334786 

0.001521 

0.123113 


Table 3.5: Contrast measurements in tissue parameter maps and weighted images between 
different brain tissues 


The first column of Figure 3.8(a) shows Ti, T 2 , and proton density weighed 
spin echo images and second column shows the computed tissue parameter maps 
using these weighted images and have better contrast than the weighted images. First 
column of Figure 3.8(b) shows weighted combinations of weighted images with 
weights calculated using Ti, T 2 , and p maps respectively. These have got an excellent 
contrast compared to weighted images similar to parameter maps. Second column of 
Figure 3.8(b) demonstrates the segmented images using weighted images, parameter 
maps, and weighted combination images respectively. Segmentation is done using a 
variant of maximum likelihood estimation (R.K.S. Rathore, S. Datta, R.K. Gupta, S.B. 
Rao and R. Verma (2001)). 


































PD weighted 


PD Map 


Figure 3.7; Histogram plots of white matter values in weighted images and in computed 
parameter maps 














PD weighted 

Figure 3.8(a): 1'' column: weighted spin echo images 2 ^^ column: corresponding computed 
parameter maps 


Ti weighted 


Ti Map 


T2 weighted 


T2 Map 











Weighted combination with PD Segmented using weighted combinations 


Figure 3.8(b): T' column: weighted combination images 2”“^ column: segmented with different 
images 









Chapter 4 

Quantitation of Cystic Cavity, Perilesional 
Gliosis, and Edema 


In this chapter three pathologies; cystic cavity, gliosis and edema, are quantitatively 
analyzed. Magnetization transfer spin-echo Ti (MT SE Ti), Ti, Ta and PD weighted 
images along with the corresponding MTR, Ti, Ta and PD computed maps are used for 
the present quantitative study. Ti, Ta, and PD maps are computed using the method 
described in subsection 2.2.4 of chapter 2. MTR values are computed as mentioned in 
subsection 1.2.4 of chapter 1. All the weighted images are normalized to have values 
from 0 to 255 so that the computed numericals on weighted images could be compared. 
Parameter values are used without any normalization because of their significance and 
existing standards of quantitative MRI in medicine. Even though it is customary to 
represent PD computed values with respect to CSF, in the present study these values are 
presented as they are obtained. First, each pathological region is segmented using the 
multi-resolution MLE algorithm (R.K.S. Rathore, S. Datta, R.K. Gupta, S.B. Rao, and 
R. Verma (2001)) with Ti, T 2 , and PD parameter maps to obtain the required region of 
interest (ROI). For comparison purpose, the contra-lateral normal region is selected by 
carefully placing the same size and shape of the segmented pathological region on the 
opposite hemisphere of the brain on each slice. This is done to have almost same 
number and type of tissue components in the pathological region as well as in the 
normal region. Otherwise, the comparison would not be justified. In some cases it may 
not be possible to select in this manner. For example, in the case of edema, in one of the 
slices, the pathological region has pushed the lateral ventricle and if one chooses the 
contra-lateral normal region in same location on opposite hemisphere of the brain, it 
would contain sufficient amount of contribution from CSF in the ventricle and 
consequently would contaminate the quantitative values. In such critical cases, keeping 
physiological structures in view, a small part of the segmented pathological region is 
ignored to obtain more meaningful contra-lateral normal region. 



For the present study, MR imaging of the brain is performed on 1.5T 
superconducting system using circularly polarized head coil. Conventional spin-echo 
proton density, T 2 (TR/TE/l,2/n=2200/12,80/l), and Ti (1012/14/2) weighted imaging 
is performed in axial plane using 256x256 matrix size, 0.1 mm inter-slice gap and 5mm 
slice thickness. MT SE Ti weighted MR imaging is performed with exactly the same 
parameters as for Ti weighted except for the off-resonance pulse. 

In section 4.1, the statistical tools used for the quantitative analysis are 
discussed. Sections 4.2, 4.3, and 4.4 are devoted to the quantitative study of the 
pathologies: cystic cavity, perilesional gliosis, and edema with respect to the contra- 
lateral normal regions. 

4.1. Statistical Tools 

Following W.H. Press, S.A. Tenkolsky, W.T. Vellerling, and B.D. Flannery (1998), in 
subsections 4.1.1 and 4.1.2, F-test and Student's t test and their probability functions are 
briefly described. These are used to study the significance of differences in variances 
and means in pathological regions compared with contra-lateral normal regions. In 
subsection 4.1.3, formula for computing linear correlation coefficients is given. 

4.1.1. F-test for Significantly Different Variances 

The F-test tests the hypothesis that two samples have different variances by trying to 
reject the null hypothesis that their variances are actually consistent. The statistic F is 
the ratio of one variance to the other, so values either »1 or «1 will indicate very 
significant differences. 

F-distribution Probability Function 

A certain statistic F, essentially the ratio of the observed dispersion of the first sample to 
that of the second one, is calculated. The probability that F would be as large as it is if 
the first sample's underlying distribution actually has smaller variance than the second's 
is denoted QiF\Vi,V 2 ) where Vi, V 2 are the number of degrees of freedom in the first 
and second samples respectively. In other words, QiF | v'j,V 2 ) is the significance level 
at which the hypothesis "1 has smaller variance than 2" can be rejected. A small 
numerical value implies a very significant rejection, in turn implying high confidence in 
the hypothesis " 1 has variance greater or equal 2." 



^(FlVpVj) has the limiting values 

<2(0|v,,V2) = 1 and 2(°°|v,,V2) = 0 . 
Its relation to the incomplete beta function is 

Q{F\v„v^) = l 

V2 + V1F 


4.1.2. Student’s t-test for Significantly Different Means 


Applying the concept of standard error, the conventional statistic for measuring the 
significance of a difference of means is termed Student's t. When the two distributions 
are thought to have the same variance, but possibly different means, then Student's t is 
computed as follows; First estimate the standard error of the difference of the means, 
Sd, from the "pooled variance" by the formula 



^ 1 n 

V N,+N,-2 

1 

^A ^B 

\ A B J 




where each sum is over the points in one sample, the first or second, each mean likewise 
refers to one sample or the other, and Na and Nb are the numbers of points in the first 
and second samples, respectively. 

Now the Student's t statistic can be computed from 


t = 


^A-^B 


and degrees of freedom can be computed from 

df = N,+N,- 2 . 

If two distributions have significantly different variances, the t statistic can be computed 
from 


t = - 


Xa-^b 


[Yarix,)/N,+Yax(x,)/N,r’ 
and degrees of freedom can be computed from 


df = 


Var(x^)^ Var(x^ 


iV. 


iV. 


\Yax(xJ/Nj ^ [Var(x,)/iV,r ' 


N,-l 


N,-l 


To find out whether the two data sets have variances that are significantly different, the 
F-test described in the previous section can be used. 

Student's Distribution Probability Function 


A(r I v) is the probability, for v degrees of freedom, that a certain statistic i (measuring 

the observed difference of means) would be smaller than the observed value if the 
means were in fact the same. Two means are significantly different if, e.g., 
A(r|v)>0.99. In other words, 1-A(rlv)is the significance level at which the 
hypothesis that the means are equal is disproved. 

The mathematical definition of the function is 


Mt 1 V) = f (1 + dx . 

Limiting values are 

A(0|v) = 0, A(oo|v) = l. 

A(r I v) is related to the incomplete beta function (a, b) by 

A(r|v) = l-I_^(f,l). 


4.1.3. Linear Correlation Coefficients 


Linear correlation coefficient, r for pair of quantities (x;,y,), i = 1,...,N, can be 
computed by the formula 

where x is the mean of x,- 's and y is the mean of y,- 's. 


4.2. Quantitation of Cystic Cavity 

In this section we have considered a case of young man of 35 years old with trauma in 
frontal lobe in 1997. For two years patient was asymtomatic but, after two years patient 
presented with intracable complex partial seizures with headache. MRI was done in 



of encephalomalesia and cystic cavity. This is hyperintense on T 2 and hypointense on Ti 
weighted images. On MT SE Tj weighted image, there is gliosis which is hyperintense 
and this abnormality is beyond the T 2 abnormality. Five consecutive cross-sections of 
each 5mm thickness and 0.1mm inter-slice gap are considered and these are through the 
level of body lateral ventricles to third ventricles. 

Table 4.1 shows relatively low mean values in the cystic cavity regions 
compared with the corresponding contra-lateral normal regions on MT SE Ti weighted 
images. Variances in cystic cavity are relatively low. Significantly low mean and 
variances of MTR values in cystic cavity regions could be observed in Table 4.2. 


MT SE 
Ti wtd. 

Mean ± S.D. 



t 

P 

F 

1101111111 

Normal 

Cystic cavity 

Slice 1 

75.00±5.28 



1.80e-23 

2.189308 

1.50e-24 

Slice 2 

74.1516.72 

63.8114.49 

38.567856 

O.OOe+00 

2.233288 

5.87e-33 

Slice 3 

69.2516.71 

61.6813.50 





Slice 4 


57.3914.86 

59.287617 

O.OOe+00 



Slice 5 

69.6116.99 

64.0115.77 

17.028128 

O.OOe+00 

1.466500 

1.46e-07 


Table 4.1: Mean and standard deviation of normal and cystic cavity regions and Student’s t, F 
statistics and their probabilities on MTSE Ti weighted images 


MTR 

Mean 1 S.D. 

t 

P 

F 

P 

Normal 

Cystic cavity 

Slice 1 

29.5416.96 

6.7113.35 

78.041046 

O.OOe+00 

4.324141 

O.OOe+00 

Slice 2 

29.2618.24 





O.OOe+00 

Slice 3 

25.6818.52 

7.5113.69 

50.448616 

O.OOe+00 

5.332147 

O.OOe+OO 

Slice 4 

28.4518.80 

6.5013.94 

67.571609 

O.OOe+00 

4.985687 

O.OOe+00 

Slice 5 

27.3719.34 



O.OOe+00 

3.983310 

O.OOe+00 


Table 4.2: Mean and standard deviation of normal and cystic cavity regions and Student’s t, F 
statistics and their probabilities on MTR maps 


Table 4.3 shows low mean pixel intensities in the cystic cavity regions with 
respect to the corresponding normal regions on Ti weighted spin echo images. The 
variances in cystic cavity are significantly low. In contrast to the weighted images Ti 
maps show significantly high Ti relaxation times in pathological region compared with 
the corresponding normal regions as shown in Table 4.4. 
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Ti wtd. 

Mean ± S.D. 

t 

P 

F 


Normal 

Cystic cavity 


99.29112.56 

71.7714.18 

54.964119 




Slice 2 

104.25116.25 

67.0215.14 

65.873589 

O.OOe+00 

O.OOe+00 

O.OOe+00 

Slice 3 

89.91115.28 

63.2813.86 

43.552238 

O.OOe+00 

O.OOe+00 

O.OOe+00 



57.1015.40 

76.549767 

O.OOe+00 

O.OOe+00 

O.OOe+00 

Slice 5 

96.64115.50 

71.6017.68 

39.895119 

O.OOe+00 

O.OOe+00 

O.OOe+00 


Table 4.3: Mean and standard deviation of normal and cystic cavity regions and Student’s t, F 
statistics and their probabilities on Ti weighted images 


Ti 

values 

Mean 1 S.D. 

t 

P 

F 

P 

Normal 



Cystic cavity 

Slice 1 


1970.731224.24 

-43.25634 


4.73813 

O.OOe+00 

Slice 2 

1101.011796.97 

2297.301337.48 

-41.69613 

O.OOe+00 



Slice 3 

1305.801775.07 

2095.731309.57 





Slice 4 

1017.911473.09 

2694.971424.75 

-78.36119 

O.OOe+00 

1.24057 

1.38e-03 

Slice 5 

1182.331627.67 

2069.661344.29 

-34.15867 

O.OOe-i-00 

3.32356 

O.OOe+00 


Table 4.4: Mean and standard deviation of normal and cystic cavity regions and Student’s t, F 
statistics and their probabilities on Ti maps 


Tables 4.5 and 4.6 show high mean values with less variances in the cystic 
cavity region compared with normal regions on both T 2 weighted image and T 2 map. 
Similar to Ti weighted images and Ti maps, cystic cavity has relatively low pixels 
intensities on PD weighted images and high mean values on PD maps compared with 
normal regions as shown in Tables 4.7 and 4.8 respectively. 


T 2 wtd. 

Mean 1 S.D. 

t 

P 

F 

P 

Normal 

Cystic cavity 

Slice 1 

145.21126.65 

223.87112.93 

-70.190910 



O.OOe+00 

Slice 2 

130.75126.06 

192.90114.69 

-62.642597 



O.OOe+00 

Slices 

124.50125.65 

177.3619.46 

-49.847652 

O.OOe+00 

7.343417 

O.OOe+00 

Slice 4 

115.97124.97 

167.90114.25 

-53.653759 

O.OOe+00 

3.070495 

O.OOe+00 

Slice 5 

123.69124.96 

187.48116.69 

-58.517342 

O.OOe+00^ 

2.235762 

7.66e-28 

Table 4.5: 

Mean and stand 

lard deviation of normal and cystic cavity regions and Student’s t, F 


statistics and their probabilities on T 2 weighted images 































































Q 

Mean ± S.D. 

t 

P 

F 

P 

Normal 

Cystic cavity 


72.88±24.97 

145.22116.43 






69.84±30.63 

127.20117.97 

-48.714535 





80.65±44.52 

150.82119.60 

-37.199162 


5.157619 



74.76±62.09 

131.23118.16 

-25.940105 

O.OOe-t-OO 

11.678875 

O.OOe-t-OO 


73.47±26.68 

135.65120.54 

-50.875862 


1.686646 

8.13e-13 


Table 4.6: Mean and standard deviation of normal and cystic cavity regions and Student’s t, F 
statistics and their probabilities on T 2 maps 


PD wtd. 

Mean 1 S.D, 

t 

P 

F 

P 

Normal 

Cystic cavity 


130.4716.54 

121.5514.94 

28.747776 

O.OOe-^OO 

1.757371 

1.34e-13 



129.9115.95 

41.817730 

0.00e-+00 

2.132167 

1.56e-29 


134.0819.29 

118.5013.73 

40.109089 

O.OOe-hOO 

6.176240 

O.OOe-hOO 


139.7818.95 

124.9917.44 

37.722557 

O.OOe+00 

1.444870 

5.05e-08 



133.2418.64 

14.249944 

2.45e-43 

1.572978 

5.24e-10 


Table 4,7: Mean and standard deviation of normal and cystic cavity regions and Student’s t, F 
statistics and their probabilities on PD weighted images 


PD 

values 

Mean 1 S.D. 

t 

P 

F 

P 

Normal 

Cystic cavity 

Slice 1 

1558.151107.13 

1564.29196.33 

-1.12700 

2.60e-01 

1.23685 

5.01e-03 

Slice 2 

1620.921184.10 

1665.551118.54 

-6.14876 

9.57e-10 





1600.321155.24 

-12.38100 

2.12e-33 

2.07856 

1.07e-20 



1773.321143.89 

-30.08028 



8.98e-ll 




-6.39808 

2.09e-10 

1.53924 

3.28e-09 


Table 4.8: Mean and standard deviation of normal and cystic cavity regions and Student’s t, F 
statistics and their probabilities on PD maps 


Linear correlation coefficients between different weighted images in the cystic 
cavity region and contra-lateral normal region are tabulated in Tables 4.9 and 4.10 
respectively. MT weighted values are strongly positively correlated with Ti weighted 
and PD weighted values while negatively correlated with T 2 weighted values in the 
normal regions. In the cystic cavity region, MT weighted values showed same kind of 
















































correlation with Ti, and PD weighted values as in the case of normal regions. But, with 
T 2 weighted values MT SE Ti weighted values have shown negative correlation 
compared with contra-lateral normal regions unlike in pathological regions. In the 
normal regions, T 2 weighted values are strongly negatively correlated with Tj and PD 
weighted values while T i and PD weighted values are positively correlated. In contrast 
to the normal regions, T 2 weighted values are positively correlated with Ti and PD 
weighted values in the cystic cavity region. But the Ti and PD weighted values 
maintained same correlation as in the case of normal regions. 

Tables 4.11 and 4.12 show the linear correlation coefficients between different 
tissue parameters in the normal and cystic cavity regions respectively. MTR map has 
shown negative correlations with Ti, T 2 , and, PD maps in both the regions. Ti, T 2 , and 
PD maps are positively correlated among themselves in the normal regions. In the cystic 
cavity region, Ti with T 2 and PD maps has shown positive correlation while T 2 has 
negative correlation with PD map which is not very significant. 

Graphs 4. 1-4.4 show the mean values in cystic cavity and contra-lateral normal 
regions on weighted images in consecutive five slices. On MT SE Ti, Ti, and PD 
weighted images, the mean values in cystic cavity region are less than the mean values 
in normal regions. But, on T 2 weighted image mean values in pathological regions are 
higher than in normal regions. Graphs 4.5-4. 8 show the mean values in pathological 
regions and normal regions on computed parameter maps. On Ti, T 2 , and PD maps, 
mean values in cystic cavity regions are higher compared with normal regions. But, on 
MTR map normal region showed high values than in cystic cavity region. 




















































































































































































































































































1 2 3 4 5 


Slice Number 

Graph 4.1: Mean values in cystic cavity and contra-lateral normal regions on MT SE Ti 
weighted images 



Slice Number 

Graph 4.2: Mean values in cystic cavity and contra-lateral normal regions on Ti weighted 
images 



Slice Number 

Graph 4.3: Mean values in Cystic cavity and contra-lateral normal regions on T 2 weighted 



Slice Number 


Graph 4.4: Mean values in Cystic cavity and contra-lateral normal regions on PD weighted 


imSlCTAC 




Graph 4.5: Mean values in cystic cavity and contra-lateral normal regions on MTR maps 



Graph 4.6: Mean values in cystic cavity and contra-lateral normal regions on Ti maps 



Graph 4.7: Mean values in Cystic cavity and contra-lateral normal regions on T 2 maps 



Graph 4.8: Mean values in Cystic cavity and contra-lateral normal regions on PD maps 
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4.3. Quantitation of Perilesional Gliosis 

Neurocysticercosis is the most common parasitic disease of the central nervous system 
in the developing countries with seizures as the most common manifestation. Lesion 
localization is usually done with computed tomography (CT) and/or magnetic resonance 
imaging. Recently, hyperintense signal presumably due to gliosis has been observed on 
MT SE Ti weighted image around the healed/calcified cysticercus granuloma and is 
shown to be related to epileptogenisis and poor long-term seizure control (R.K. Gupta, 
M.K. Kathuria, S. Pradhan (1999)). For the present study the same patient whose case 
was discussed in the previous section is considered. A perilesional gliosis around the 
cystic cavity is observed. The presumed perilesional gliosis is considered present when 
hyper intensity is seen on MT SE Ti weighted image around the lesion; not/barely 
visible abnormality on T 2 weighted images. In this case also we have considered five 
consecutive slices each of thickness 5mm with inter-slice gap 0.1mm which are from 
the level of body lateral ventricles to third ventricles. Even though the gliosis is slightly 
visible on PD weighted images, it is difficult to observe it on Ti weighted images. In the 
case of T 2 weighted images it has merged with the cystic cavity's hyper intensity. 

Table 4.13 shows the mean and standard deviation values of pixels in 
perilesional gliosis region and corresponding contra-lateral normal regions on MT SE Ti 
weighted images of five consecutive slices. Student's t and F statistics and their 
corresponding probabilities are also tabulated to study the significance of the changes in 
mean and variances of pixel values in gliosis region compared to the contra-lateral 
normal region. In all the slices the mean values are significantly high (very small 
probabilities of t-statistic) compared to normal region and justifies the appearance of 
gliosis on MT SE Ti weighted images as hyper-intense. The high probabilities of F- 
statistic shows that the variations with in the gliosis region are not significant. 


MT SE 
Ti wtd. 

Mean ± S.D. 

t 

P 

F 

P 

Normal 

Gliosis 

Slice 1 

75.57±1.98 

85.05±2.44 

66.625107 

O.OOe+00 

1.521773 

3.78e-06 

Slice 2 

74.82±2.34 

82.71±2.32 

48.645851 

O.OOe+00 

1.016214 

8.70e-01 

Slice 3 

71.22±2.92 

80.83±3.01 

45.799877 

O.OOe+00 

1.077150 

4.46e-01 

Slice 4 

73.03±2.96 

79.64±2.72 

24.095078 

O.OOe+OO 

1.184203 

2.16e-01 

Slice 5 

70.8G±2.50 

80.74±2.63 

50.604851 

O.OOe+00 

1.112922 

3.22e-01 
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from Table 4.16. This increment in Tj relaxation times is around 100 ms compared to 
normal region. The appearance of gliosis on Ti maps would be bright. 


Ti 

values 

Mean ± S.D. 

t 


F 

P 

Normal 

Gliosis 

Slice 1 

822.891180.06 






Slice 2 

800.991235.30 

993.221167.86 

13.548338 

6.33e-38 

1.964858 

9.64e-12 

Slice 3 

859.911200.25 

942.641163.82 

6.761215 

2.55e-ll 

1.592198 

1.82e-06 

Slice 4 






9.09e-19 

Slice 5 

916.541217.73 

1041.311204.7 

7.731840 

3.78e-14 

1.131243 

2.54e-01 


Table 4.16: Mean and standard deviation of normal and gliosis regions and Student’s t, F 
statistics and their probabilities on Ti maps 


Tables 4.17 and 4.18 show significantly high values in gliosis compared to 
normal region on T 2 weighted and T 2 maps respectively. On both T 2 weighted and T 2 
maps gliosis would appear as bright. In both the cases the variations in gliosis region are 
significant. 



Mean 1 S.D. 

T 

P 

F 

P 



Slice 1 

126.44113.15 



O.OOe+00 

6.016067 

O.OOe+00 

Slice 2 

118.82113.65 

168.52126.63 

33.821072 


3.804779 

1.89e-39 

Slice 3 

108.60116.80 

148.95127.20 

25.388273 



1.52e-17 

Slice 4 

107.94120.49 

140.44119.41 

16.879768 

O.OOe+00 

1.113660 

4.31e-01 









Table 4.17: Mean and standard deviation of normal and gliosis regions and Student’s t, F 
statistics and their probabilities on T 2 weighted images 


T 2 

values 

Mean 1 S.D. 

T 

P 

F 

P 

Normal 

Gliosis 

Slice 1 

62.6215.34 

86.61118.45 

27.616093 

O.OOe+00 

11.938922 

O.OOe+00 

Slice 2 

61.5216.23 

81.24113.52 

26.987143 

O.OOe+00 

4.710935 

O.OOe+00 

Slice 3 

63.9018.52 

83.25115.86 

21.706852 

O.OOe+00 

3.098142 

1.97e-29 

Slice 4 

63.07110.57 

76.5319.31 

13.997484 

5.99e-37 

1.288386 

6.39e-02 

Slice 5 

62.8717.49 

85.64115.02 

25.121229 

O.OOe+00 

4.017428 

• j 1 

2.62e-35 

n.*. a. 


Table 4.18: Mean and standard deviation of normal and gliosis regions and Student’s t, F 







































Table 4.19 shows significantly high values in gliosis region on PD weighted 
images confirming its hyper intensity on these weighted images. The differences in 
variations between gliosis and normal regions are not significant. Table 4.20 also shows 
significant increment in computed PD values in gliosis region with no significant 
difference in variances. Ta weighted image values showed more variation compared to 
Ti, PD, and MT weighted image values within the gliosis region. 


PD 

wtd. 

Mean ± S.D. 

t 

P 

F 

P 

Normal 



126.9613.68 

134.4214.16 

29.673491 


1.273812 

7.56e-03 

Slice 2 

118.82113.65 

153.9514.16 



1.179972 

9.22e-02 

Slice 3 

108.60116.80 

144.4215.40 

28.110928 

O.OOe+00 

1.006975 

9.42e-01 

Slice 4 


150.7515.46 

17.827040 


1.369388 

2.16e-02 







9.54e-01 


Table 4.19: Mean and standard deviation of normal and gliosis regions and Student’s t, F 


statistics and their probabilities on PD weighted images 


PD 

values 

Mean 1 S.D. 

t 

P 

F 

P 

Normal 

Gliosis 

Slice 1 

1472.19171.80 

1518.44160.61 

10.885725 

3.96e-26 

1.403464 

1.86e-04 

Slice 2 

1516.30182.16 

1628.12161.28 

22.225010 

O.OOe+00 

1.797602 

3.07e-09 

Slice 3 

1508.73181.85 

1581.35168.45 

13.816566 

2.63e-39 

1.438682 

1.86e-04 




10.582300 

2.05e-23 

1.993117 

5.68e-07 



1627.31166.94 

11.161126 

1.07e-26 

1.564804 

3.69e-05 


Table 4.20: Mean and standard deviation of normal and gliosis regions and Student’s t, F 
statistics and their probabilities on PD maps 


Linear correlation coefficients between different weighted images in normal and 
gliosis regions are shown in Tables 4.21 and 4.22 respectively. Normal regions on MT 
SE Ti weighted images did not show consistent correlations with Ti and T 2 weighted 
images. However, it has positive correlation with PD weighted images. Ti weighted 
images in the normal regions have shown negative correlation with T 2 and PD weighted 
images while T 2 and PD are positively correlated. In gliosis region, MT SE Ti weighted 
images are mostly positively correlated with T 2 and PD weighted images and with Ti 











































weighted images are negatively correlated with T 2 and PD weighted images in gliosis 
regions. T 2 and PD also have shown similar kind of correlation in gliosis regions 
compared with normal regions. 

Table 4.23 shows linear correlation coefficients between MTR, Tj, T 2 and PD 
maps in the normal regions. MTR is significantly negatively correlated with Ti, T 2 and 
PD values. MTR is more strongly correlated with Ti map than with T 2 and PD maps. Ti, 
T 2 , and PD are all significantly positively correlated among themselves. The correlation 
between T 2 and PD is relatively less compared with the correlations between Ti , T 2 and 
T,,PD. 

Table 4.24 shows linear correlation coefficients between different computed 
parameter maps in gliosis regions. All correlations are similar to the normal regions. 

Graphs 4.9-4.12 show the mean values in perilesional gliosis and contra-lateral 
normal regions on weighted images in consecutive five slices. On MT SE Ti, T 2 , and 
PD weighted images, the mean values in gliosis regions are relatively high compared 
with normal regions as shown in Graphs 4.9, 4.11, and 4.12. Graph 4.10 shows no 
significant difference in mean values between pathological and normal regions on Tj 
weighted images. Graph 4.13-4.16 show the mean values in gliosis and normal regions 
on computed parameter maps. Graph 4.15 and 4.16 show, respectively, the increase in 
mean values on T 2 and PD maps in gliosis compared with normal regions while Graph 
4.13 shows decrease in mean MTR values. In contrast to Ti weighted mean values. 
Graph 4.14 shows a significant increase in mean computed Ti values in pathological 
regions. 

This analysis indicates of the usefulness of magnetization transfer effect for the 
diagnosis of perilesional gliosis. The observation of significant high Ti values in gliosis 
region could be used for the applicability of other pulse sequences which enhance the 
contrast of this region for diagnosis purposes. For example, in section 3.5, synthesized 
IR images at parameters TR/TE/TI = 6000/40/2500 ms have shown an excellent 
contrast between the gliosis and surrounding brain tissue. 
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Graph 4.9: Mean values in perilesional gliosis and contra-lateral normal regions on MT SE 
Ti weighted images 



Graph 4.10: Mean values in perilesional gliosis and contra-lateral normal regions on Ti 
weighted images 



Graph 4.11: Mean values in perilesional gliosis and contra-lateral normal regions on T2 
weighted images 



weighted images 



Mean T. values 



Graph 4.13: Mean values in perilesional gliosis and contra-lateral normal regions on MTR maps 



Graph 4.14: Mean values in perilesional gliosis and contra-lateral normal regions on Ti maps 




Graph 4.16: Mean values in perilesional gliosis and contra-lateral normal regions on PD maps 


4.4. Quantitation of Edema 


In the present section a quantitative analysis of edema has been done in the case of a 
single patient. This patient is a 45 years old female who was in good health till April 
1999 and was diagnosed with a tapeworm infection in July 1999. In July 1999, she had 
experienced sudden onset of seizures which manifested initially as difficulty with 
speech and secondary generalization, and later with left side headache and with simple 
partial seizures. MRI was done in July 1999 and there was healing lesion in left parietal 
lobe with extensive edema and anti epileptic drugs (AED) were started. For the present 
study, eight consecutive slices each of them with 5mm slice thickness and 0.1mm inter- 
slice gap are considered. These cross-sections are from high parietal lobe to the level of 
third ventricles. There is single lesion in left high parietal region with extensive edema 
which is extending from high parietal lobe to parieto-occipital lobe. 

Table 4.25 shows a significant increment in MT weighted values in edema 
regions in comparison with contra-lateral normal regions. But there is no significant 
difference in variances in both the regions. In contrast to this, MTR values are 
significantly low in edema regions except in slice 8 as shown in Table 4.26. 


MT SE 
Ti wtd. 

Mean ± S.D. 

t 

P 

F 

P 

Normal 

Edema 

Slice 1 

64.1414.44 

68.7113.53 

-22.338718 

O.OOe-i-00 

1.580010 

2.52e-10 

Slice 2 

66.9114.97 

72.2313.10 

-27.825333 

O.OOe+00 

2.562336 

1.40e-45 

Slice 3 

70.0313.73 

75.0412.80 


O.OOe-t-OO 

1.778403 

3.34e-27 

Slice 4 


75.8712.09 

-59.428463 

O.OOe+00 

2.004717 

1.81e-35 

Slice 5 

72.0312.22 

77.5312.38 

-64.569191 

O.OOe-t-00 

1.147881 

8.27e-03 

Slice 6 

71.5712.91 

74.7411.54 


O.OOe+00 

3.534718 

O.OOe+00 

Slice 7 

67.2612.71 

76.2911.63 

-71.200199 

O.OOe-i-00 

2.764753 

2.36e-35 

Slice 8 



-38.699139 

O.OOe-tOO 

4.359744 

O.OOe+00 


Table 4.25: Mean and standard deviation of normal and edema regions and Student’s t, F 
statistics and their probabilities on MTSE Ti weighted images 





















MTR 

Mean ± S.D. 

t 

P 

F 

P 

Normal 

Edema 

Slice 1 

28.44±7.53 

21.4414.49 





Slice 2 

35.82±6.25 


33.725716 




Slice 3 

27.55±6.87 

18.1214.15 

44.336849 




Slice 4 

37.59±4.88 

27.9712.60 

62.823109 


3.515889 

O.OOe+00 

Slice 5 

29.0415.90 

16.4813.11 

72.119064 



O.OOe+00 



28.4211.62 

84.767334 

o.ooe+00 

3.505460 

O.OOe+00 

Slice 7 

29.1814.36 

17.1313.42 



1.622852 

1.96e-09 

Slice 8 

36.2917.33 

35.8812.93 

1.353652 

1.76e-01 

6.267653 

O.OOe+00 


Table 4.26: Mean and standard deviation of normal and edema regions and Student’s t, F 
statistics and their probabilities on MTR maps 


On Ti weighted image, edema has significantly low values and Ti map shows 
high values as shown in Tables 4.27 and 4.28 respectively. Both on T 2 weighted images 
and T 2 maps edema, on overall, has significantly high values compared with contra- 
lateral normal regions with significant difference in variances as demonstrated in Tables 
4.29 and 4.30 respectively. 

From Tables 4.3 1 and 4.32, on overall, PD weighted images and PD maps show 
significantly high values in edema region compared to the contra-lateral normal regions. 


T, 

wtd. 

Mean 1 S.D. 

t 

P 

F 

P 

Normal 

Edema 

Slice 1 

90.57110.69 

87.6515.04 

6.863605 

9.69e-12 

4.496512 

O.OOe+00 

Slice 2 

97.34111.25 

92.9714.46 

11.083804 

1.08e-27 

6.371774 

O.OOe+00 

Slice 3 

96.7519.80 

91.1114.01 

20.118578 

O.OOe+00 

5.965241 

O.OOe+00 

Slice 4 

103.7918.43 

97.1814.02 

25.547007 

O.OOe+00 

4.393178 

O.OOe+00 

Slice 5 

101.8918.25 

92.7114.45 

37.519688 

O.OOe+00 

3.427851 

O.OOe+OO 

Slice 6 

111.0115.70 

98.3011.83 

57.475475 

O.OOe+00 

9.607631 

O.OOe+00 

Slice 7 

95.4516.33 

92.3513.88 

10.393809 

2.54e-24 

2.657025 

7.22e-33 

Slice 8 

109.1712.81 

99.82113.80 

17.47448 

O.OOe+00 

24.070436 

• . 1 

O.OOe+OO 

C \ j . 1 4.5 « + 


Table 4.27: Mean and standard deviation of normal and edema regions and Student s t, F 


statistics and their probabilities on T i weighted images 






























































T, 

values 

Mean ± S.D. 

t 

P 

F 

P 

Normal 

Edema 

Slice 1 

969.72±403.45 

1202.631158.48 

-14.93623 

O.OOe-^00 

6.48048 

O.OOe+00 

Slice 2 

741.26±237.30 

1053.141109.95 

-36.63475 

O.OOe+00 

4.65783 

O.OOe+00 

Slice 3 

950.40±343.78 

1336.091132.85 

-39.59203 

O.OOe+00 

6.69596 

O.OOe+OO 

Slice 4 

719.221197.55 

1065.89189.40 

-57.73639 

O.OOe-hOO 

4.88250 

O.OOe+OO 

Slice 5 

921.671247.87 

1395.361141.77 

-63.57553 

O.OOe+00 

3.05685 

O.OOe+OO 

Slice 6 

610.531119.80 

1082.77169.07 

-92.49870 

O.OOe+00 

3.00832 

O.OOe+OO 

Slice 7 

904.801178.70 

1248.231102.53 

-41.53815 

O.OOe+00 

3.03800 

1.42e-41 

Slice 8 

776.31189.17 

828.461473.23 

-2.85256 

4.40e-03 

28.16275 

O.OOe+OO 


Table 4.28: Mean and standard deviation of normal and edema regions and Student’s t, F 


statistics and their probabilities on Ti maps 



Mean 1 S.D. 

t 

P 

F 

P 

1 2 Wlu. 

Normal 

Edema 

Slice 1 

145.04120.12 

215.10125.32 

-60.174534 

O.OOe+OO 

1.583861 

2.04e-10 

Slice 2 

121.75113.28 

198.09113.26 

-124.085564 

O.OOe+OO 

1.002524 

9.69e-01 

Slice 3 

135.02119.93 

220.51113.75 

-133.400696 

O.OOe+OO 

2.100069 

1.06e-43 

Slice 4 

124.87113.53 

215.2517.73 

-209.274643 

O.OOe+OO 

3.060679 

O.OOe+OO 

Slice 5 

123.99114.52 

209.5717.61 

-200.056305 

O.OOe+OO 

3.637607 

O.OOe+OO 

Slice 6 

120.59110.77 

211.6717.70 

-186.219406 

O.OOe+OO 

1.956700 

2.17e-19 

Slice 7 

129.65112.26 

215.6017.10 

-151.130096 

O.OOe+OO 

2.981824 

2.64e-40 

Slice 8 

115.33116.58 

166.08121.05 

-49.834938 

O.OOe+OO 

1.610588 

4.35e-10 

Table 4.29: Mean and standard deviation of normal and edema regions and Student’s t, F 
statistics and their probabilities on T 2 weighted images 

n 

Mean 1 S.D. 

t 

P 

F 

P 

Normal 

Edema 

Slice 1 

68.39124.77 

95.75115.86 

-25.839197 

O.OOe+OO 

2.439183 

3.55e-34 

Slice 2 

67.83121.64 

106.47113.15 

-46.783211 

O.OOe+OO 

2.707615 

O.OOe+OO 

Slice 3 

70.22125.60 

111.56111.78 

-55.477112 

O.OOe+OO 

4.723387 

O.OOe+OO 

Slice 4 

65.88111.70 

116.8018.50 

-127.08400 

O.OOe+OO 

1.896296 

2.28e-30 

Slice 5 

67.46112.06 

120.0419.29 

-132.36201 

O.OOe+OO 

1.685795 

2.55e-23 

Slice 6 

63.9914.68 

116.4916.04 

-185.91993 

O.OOe+OO 

1.663032 

7.46e-12 

Slice 7 

68.7918.72 

118.0218.98 

-97.939568 

O.OOe+OO 

1.059855 

4.69e-01 

Slice 8 

67.49120.03 

90.04113.06 

-24.825291 

O.OOe+OO 

2.352493 

1.35e-28 


Table 430: Mean and standard deviation of normal and edema regions and Student s t, F 















































Mean ± S.D. 


t 


PD 


wtd. 

Normal 

Edema 

t 

F 


P 

Slice 1 

128.40±10.68 

137.0913.27 







133.9614.06 

-30.867407 


6.996948 


Slice 3 

129.05±9.63 


-41.299812 




Slice 4 

127.87±5.84 

138.6513.54 



2.722035 

O.OOe+00 

Slice 5 

135.8415.34 

145.6613.85 

-57.155006 

o.ooe+00 

1.923388 

1.99e-35 

Slice 6 





1.534128 

7.84e-09 

Slice 7 

114.2716.12 

124.9212.55 

-39.985962 


5.766409 




147.2613.49 

-41.327747 

O.OOe-fOO 

4.439348 

O.OOe+OO 


Table 4.31: Mean and standard deviation of normal and edema regions and Student’s t, F 
statistics and their probabilities on PD weighted images 


PD 

values 

Mean ± S.D. 

t 

P 

F 

P 

Normal 

Edema 

Slice 1 

1 170.761145.66 

1222.27157.09 

-9.152752 

1.70e-19 

6.509777 


Slice 2 

1139.661124.56 

1216.29155.30 

-17.277565 

O.OOe+00 


O.OOe-i-OO 

Slice 3 



-21.853928 


8.975176 

0.00e-(-00 

Slice 4 

1214.56184.06 

1276.24142.89 

-23.603373 

O.OOe-rOO 



Slice 5 

1309.91193.98 

1408.28156.37 



2.743973 


Slice 6 

1253.86144.29 

1341.01143.05 

-38.209717 

O.OOe-i-OO 

1.058820 ^ 

4.39e-01 

Slice 7 

1334.341101.08 

1435.97140.40 

-23.264595 

0.00e-(-00 



Slice 8 

1293.591106.13 

1330.06131.78 

-8.668811 

1.20e-17 

11.153337 

O.OOe+00 


Table 4.32: Mean and standard deviation of normal and edema regions and Student’s t, F 
statistics and their probabilities on PD maps 


Tables 4.33 and 4.34 show the linear correlation coefficients among MT SE Ti, 
Ti, T 2 , and PD weighted images in the normal regions as well as in the edema regions. 
In both the regions MT SE Ti weighted values have shown positive correlation with Ti 
and PD weighted values. But with T 2 weighted images no consistency in correlations is 
observed. Except in slice 6, Ti weighted values have shown negative correlation with Ta 
weighted values in both pathological and normal regions. However, the correlation 
between Ti and T 2 weighted images in slice 6 are not significant. The correlations 

between T2 and PD weighted values have not shown consistency. 











































A similar study done on parameter maps is tabulated in Tables 4.35 and 4.36. 
MTR maps have shown negative correlation with Ti, T 2 , and PD maps in both edema 
and normal regions. With Ti map the correlations are significant than with T 2 and PD 
maps. Ti maps have shown positive correlation with T 2 and PD maps in both the regions 
while T 2 maps are negatively correlated with PD maps. 

Graphs 4.17-4.20 show the mean values in edema and contra-lateral normal 
regions on weighted images in consecutive eight slices. On MT SE Ti, and PD weighted 
images, the mean values in edema regions are relatively high compared with normal 
regions as shown in Graphs 4.17, 4.19, and 4.20. Graph 4.18 shows relatively low mean 
values in edema regions. Graphs 4.21-4.24 show the mean values in edema and normal 
regions on computed parameter maps. On Ti, T 2 , and PD maps, mean values in edema 
regions have shown relatively high values compared with normal regions as shown in 
Graphs 4.22, 4.23, and 4.24. But, Graph 4.21 shows the decrease in MTR values in 
edema regions. 

We suggest that the strong quantitative mean differences and correlations may be 
investigated by the biochemists/pathologists for the underlining cases so that the 
quantitations could be related with different stages of the pathology. As MT requires 
much more careful attention and MTR is susceptible to patient motion, quantitations 
based on routine images ought to be very useful diagnostic tools and may be 
investigated using a large data base. 
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Parametric I r(MTR.T,) I rCMTR.T,) r(MTR. PD) KT-.T^ ) !:gj^ P ) ^(^ 2, PD ) 

normal Normal P Normal P Normal P Normal P Normal P Normal P 

regions — — — u,,. , 

Slice 1 -0.7891 O.OOe+00 -0.4698 1.19e-43 -0.1477 3.78e-43 0.3120 6.77e-19 0.4789 1 .406-45 -0,5716 O.OQe+00 

Slice 2 -0.8150 O.OOe+00 -0.4514 O.OOe+00 0.1012 2.20e-03 0.1793 4.87e-Q8 0.3048 4.36e-21 -0.6600 O.OOe+00 



Slice 5 -0.6608 O.OOe+00 -0.5354 O.OOe+00 - 0.1123 1.58e-05 0.3221 /.74e-37 0.6190 O.OOe+00 -0.3702 O.OOe+00 

Slice 6 -0.6108 O.OOe+00 -0.1312 3.59e-04 -0.4067 1.19e-30 0.0139 7.07e-01 0.8551 O.OOe+00 -0.2048 2.08e-08 

Slice 7 -0.7824 0 OOe+00 -0.5762 O.OOe+00 -0.0040 9.02e-01 0.4113 8.58e-27 0.4339 5.89e.30 -0.4921 2.98e-39 

Slices -0.8485 O.OOe+00 -0.8232 O.OOe+00 -0.3403 2.64e-20 0.8110 O.OOe+00 0.6372 O.OOe+00 -0.2157 9.17e-09 

Table 4.36: Linear correlation coefficients between different parameter maps in edema regions 
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Graph 4.19: Mean values in edema and contra-lateral normal regions on T 2 weighted 
images 



Graph 4.20: Mean values in edema and contra-lateral normal regions on PD weighted 
images 
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Graph 4.21: Mean values in edema and contra-lateral normal regions on MTR map 
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Graph 4.22: Mean values in edema and contra-lateral normal regions on Ti map 
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Graph 4.23: Mean values in edema and contra-lateral normal regions on Ta map 
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Graph 4.24: Mean values in edema and contra-lateral normal regions on PD map 
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